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Summary 


This Annual Report illustrates the work carried out during the last grant-year activity on the 
Test of the Equivalence Principle in an Einstein Elevator. The activity focused on the following 
main topics: (1) analysis and conceptual design of a detector configuration suitable for the flight 
tests; (2) development of techniques for extracting a small signal from data strings with colored 
and white noise; (3) design of the mechanism that spins and releases the instrument package 
inside the cryostat; and (4) experimental activity carried out by our non-US partners (a summary 
is shown in this report). 


The analysis and conceptual design of the flight-detector (point 1 ) was focused on studying 
the response of the differential accelerometer during free fall, in the presence of errors and 
precession dynamics, for various detector’s configurations. The goal was to devise a detector 
configuration in which an Equivalence Principle violation (EPV) signal at the sensitivity 
threshold level can be successfully measured and resolved out of a much stronger dynamics- 
related noise and gravity gradients. A detailed analysis and comprehensive simulation effort led 
us to a detector’s design that can accomplish that goal successfully. 


The techniques developed for signal extraction (point 2) was focused on extracting an EPV 
signal at the threshold level from the noise affecting the detector’s output data. Techniques were 
developed for signal extraction from colored and white noises. Moreover, algorithms were 
developed/utilized for evaluating the confidence level in detecting the signal and the statistical 
error in estimating its amplitude. 


The engineers in our team worked on the design (point 3) of the mechanism that spins and 
releases the instrument package inside the cryostat during the fall of the capsule. Engineering 
drawings of the mechanism are presented with a list of key components of the mechanism and a 
description of its functions. 


Finally, a brief summary of the experimental activity (point 4) carried out during the last year 
by our non-US partners at the IFSI laboratory (Rome, Italy) is presented. The experimental work 
focused on the construction of a second differential accelerometer prototype that improves upon 
the performance of the previous prototype and its initial calibration tests in the laboratory. 
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Brief Description of Experiment 

The scientific goal of the experiment is to test the equality of gravitational and inertial mass 
(i.e., to test the Principle of Equivalence) by measuring the independence of the rate of fall of 
bodies from their compositions. The measurement is accomplished by measuring the relative 
displacement (proportional to the acceleration) of falling bodies of different materials which are 
the proof masses of a differential accelerometer 1 . 

A differential displacement at the expected frequency between the two proof masses different 
from zero will indicate a violation of the Equivalence Principle. The goal is to measure the 
Eotvos ratio 6g/g (differential acceleration/common acceleration) with a targeted accuracy of a 
few parts in 10' 15 , which is about two orders of magnitude better than the state of the art" 
Thanks to the improved accuracy the test will probe the territory where possible violations of the 
Equivalence Principle have been postulated ,v v . 

The experiment is carried out in vertical free fall with the experimental capsule released from 
a stratospheric balloon at an altitude of 40 km v '. The capsule, that is evacuated and at cryogenic 
temperature, is released from balloon and immediately afterward the detector is released inside 
the falling capsule from its top. Depending on the ballistic coefficient of the capsule, the 
detector takes between 25 s and 30 s to span the 2-m length of the capsule while the capsule falls 
between 3.1 km and 4.4 km. This technique provides an environment that is almost free of 
perturbations external to the detector. 
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Flight Detector Analysis and Conceptual Design 

Introduction 

A successful design of the differential acceleration detector for the balloon flights involves a 
careful analysis of the motion of the instrument package (that houses the detector) in free fall. 
This issue is also present in space-based detectors operating in a spacecraft that is not perfectly 
fixed in inertial space. Since rotation of the spacecraft with respect to inertial space is a highly 
desirable feature for helping in separating the violation signal from noise sources, addressing the 
issue of measuring differential accelerations in a moving detector is of general importance for 
tests of the Equivalence Principle both in free fall from balloons and in space. 

The dynamics analyses that we have conducted on various detector configurations have 
shown that the details of relative motion of the proof masses with respect to the instrument 
package (or spacecraft for a space-based test) are important for enabling the extraction of a 
violation signal from the much greater dynamics-related noise v ". 

Let us first consider a configuration in which the proof masses can only translate (i.e., they 
are constrained to move along a line) with respect to the instrument package. It is obvious that if 
the centers of mass (CMs) of the proof masses were both perfectly coincident with the CM of the 
instrument package then the dynamical coupling between the rotation of the instrument package 
and the motion of the proof masses would be null. However, this is not a realistic case because 
construction imperfections and material inhomogeneities give rise to CMs that are not 
coincident. Even if the proof masses were controlled in order to reduce the error between the 
two CMs locations, the error with respect to the CM of the instrument package or spacecraft 
would be present if the spacecraft were to have fluids on board because the spacecraft CM would 
not be perfectly stable. In other words, we must always assume that those errors are present and, 
consequently, we need to quantify their impact on our ability to extract a differential acceleration 
signal at the desired threshold level from the prevailing dynamics noise. 

In a rotational detector configuration, in which the two proof masses are constrained to rotate 
with respect to the instrument package (or spacecraft) the importance of coincidence among the 
CMs of the proof masses and the CM of the instrument package is diminished with respect to the 
importance of the equality between the moments of inertia between the two proof masses (see 
later on). 

In yet another detector configuration (that we call roto-translational), the proof masses are 
hinged off their CMs (e.g., at their edges). In this case, the motion of the proof masses with 
respect to their housing (instrument package or spacecraft) is a combination of a rotation and a 
translation and, consequently, both sets of issues pointed out previously apply to this case. 

In Annual Report #1 (Ref. v,n ) for this grant we have already analyzed the dynamics of the 
translational configuration for the flight detector. The conclusion of that analysis is that the 
precession frequency of the instrument package in free fall should be in a particular frequency 
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window in order for the dynamics-related noise to be separated from the frequency (equal to the 
rotation of the package) of the violation signal and not to saturate the differential accelerometer 
response. One problem with the translational configuration for our detector, which has 
mechanically constrained proof masses, is the difficulty in manufacturing mechanical flexural 
constraints that provide low frequency for the translational motion along the sensitive axis. 
Low-frequency flexural constraints need to be long (leading to a bulky detector) and thin 
(leading to a weak detector). We need to utilize the rotational degree of freedom for the sensing 
motion of the proof masses in order to solve this issue. Moreover, gravity gradients play a major 
role in translational configurations while their role can be strongly diminished (see later on) in 
rotational configurations of appropriate design. 

We will analyze both the roto-translational configuration and a purely rotational 
configuration of the detector in the following sections. 


Detector dynamics analysis 
Roto-translational configuration 

Motion of proof masses with respect to detector’s housing 

Let us consider a set of three rigid bodies A, B, C constrained to each other in such a way 
that the relative motion of A with respect to C is a rotation around a pivot axis <d A > and the 
relative motion of B with respect to C is a rotation around a pivot axis <dB> with <d A > and <da> 
parallel but not necessarily coincident. The pivot axes are not located at the centers of mass 
(CMs) of the proof masses but rather at distances b A and bs from the CMs of A and B. Let A 
and B be linked with C through a torsional spring dashpot system around the axes <d A > and 
<ds> characterized by elastic and damping coefficients (k A , | A ) and (ks, §b) respectively. An 
example of relationship between spin and pivot axis for one proof mass is shown in Figure 1. 


Pivot axis 



\ Spin axis 

Figure 1 Schematic of one proof mass in a roto-translational configuration. 
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Let <x,y,z> be a reference system centered at the global center of mass CM of the system 
composed by the three rigid bodies. Let y be parallel to <dA> and <ds> (such axis is exiting the 
page in Figure 2), z parallel to the central symmetry axis of body C ,and x right-handed with y 
and z. 



Since the mutual rotations between the three bodies will be of the order of a few arc seconds, 
we will assume that the components of vectors and inertia matrices with respect to the body axes 
of A and B are equal to the same components with respect to <x,y,z> (i.e., the body axes of C). 

Let cd c = (co x ,to y ,u) z ) T be the angular velocity vector of body C and consequently of the 
reference frame <x,y,z> . The angular velocity of body A and B, projected on <x,y,z> are then: 

<» A -|o) V ,0)y+^,0).j (1.1) 


0) B -|a) v ,0) v + ^-,0) z j (1.2) 


The equations of motion of the system written with respect to an inertial system centered at 
the Earth are: 


w »a(Rcm + Pa) = F 4+ f a+ f ! + F A + "LC a 


EPV 


m 


B&CM + Pb) ■ F J + F # + F f + F ! 


( 2 . 1 ) 

(2.2) 


m CPC - ~ m ApA - m BpB (2.3) 

where p is the radius vector of each individual center of mass relative to the overall CM and 
projected on the inertial system, R CM is the radius vector of CM from the earth center projected 


SAO Annual Report #2, NASA Grant N AG3-288 1 


on the inertial system, F ? is the earth gravity field force, F k the equivalent force of the internal 
spring and F ? the equivalent force of the internal damper acting on each mass. F b is the force 
transmitted to each mass from the pivotal constraint as a consequence of the rigid body attitude 
motion of A and B. We also added the acceleration a bPV acting on the mass A to simulate a 
violation of the equivalent principle. The equations of motion of the system center of mass are: 


R 


Be r 


CM 


CM 


Km 



(3) 


where F c is a perturbing force acting on the external case during free fall and m T the total mass 
of the system. 

The gravitational forces appearing in Eqs. (2. 1 ) and (2.2) can be written in terms of a gravity 
gradient matrix [T E ] with respect to CM as follows: 



(4.1) 


(4.2) 


The earth gravity gradient matrix, in inertial coordinates, has the diagonal form: 

M- 


*11 o o 
0 #22 0 
0 0 #33 


After plugging Eqs. (3) into Eqs. (2) we obtain: 

PA-[K]p A + 


F A + F f + F A F - 


m A 


^ + af v 


mj 


Pb 




Pb + 


f| + f| + f£ H 


m B 


ftlj 1 


(5) 


(6.1) 


( 6 . 2 ) 


Pc 


m A 

m c 


m B 

Pa ~Pb 

m c 


(6.3) 


The previous equations can be written with respect to the non-inertial system <x,y,z> as 
follows: 
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[ 1 F k 

r ^J r /t + — 

ffi + M r s + [o>][o)]r fi + 2[o)]r B - [r£,,]r 5 + 


+ f! + Fa 

m A 


F A epv 

— + a a 

m A 


f|+f|+f|__f £ 

m B m B 


(7.1) 

(7.2) 


*c 



(7.3) 


where r A ,r B ,r c are the radius vectors of each individual center of mass relative to CM and 
projected on the body axes <xyz> and F^, Fg are possible perturbation forces acting inside the 
detector on the proof masses. 

[R ib ], [R b ,] = rotation matrixes from the inertial to the <xyz> frame and viceversa. 

[a>] - [R/b][Rb/] = angular velocity matrix of body reference frame, 
while the gravity gradient matrix in body axes is; 


[ri-]-[R/4 r£ J*s'] 


8 xx 8xy 8xz. 
8xy 8yy 8yz 
8xz 8yz 8zz 


( 8 ) 


For defining the orientation of the body axes with respect to the inertial axes, we adopt the 
rotational Euler’s rotation sequence 2-1-3 that leads to relatively simple formulas for the external 
torques in body axes for a body spinning about an horizontal axis. Moreover, this rotation 
sequence avoids the trigonometric discontinuity for attitude oscillations about an horizontal axis. 
The physical significance of the angles in the 2-1-3 sequence with the Y inertial axis along the 
local vertical are <j> = azimuth, 0 = elevation and = spin. After defining c = cosine and s = sine, 
the 2-1-3 rotation matrix yields: 


'll 

r \2 

r n 


cl/tcf + sipsOstp slpcO 

—expsf + stpxOap 

r 2\ 

'22 

r 23 

- 

-slpap + c\psds<p cipcO 

sips <p + crpsdaj) 

431 

'*32 

*53. 


cOs<p -sO 

cOcf 


The inverse rotation (i.e., from inertial to body axes) is [/? g/ ]-[fi/gJ thanks o the 
orthogonality o rotation matrices. After introducing Eqs. (5) and (9) in Eq. (8) we obtain: 
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8 xx " r \\8\ 1 + '12822 + '13833 
8xy " r n r 2\8\ 1 + '12'22822 + ^ 3'23833 
8xz " r l l r 3l8ll + n 2 r 2 2822 + 'l3 r 33833 
gyx - r 2\ r \ l8l l + r 22 r l 2 g 22 + r 23 r l 3 g 33 

8yy " r 2l8l 1 + '22822 + r 23833 ^ ,0 ^ 

8y; ■ ^1^1*11 + '22 r 32822 + r 23 r 33833 

8yr “ '3l'l l8l 1 + '32'12822 + '33'l3833 
8jy “ r 31 r 2l8l 1 + '32 r 22822 + r 33 r 23833 

2 2 2 
8zz “ r 3l8l 1 + r 32822 + '33833 

By choosing the inertial frame such that the Y axis is parallel to R CM , the components of the 
Earth’s gravity gradient in inertial coordinates are: 


8ll 


GMe 

r CM 


8 22 


2 GMe 

r cm 


( 11 ) 


833 


GM £ 
r cm 


The forces F* and F § transmitted by the rotational spring dash-pot system can be obtained as 

shown in the following. For a displacement x A of the center of mass of A relative to C the spring 
dash-pot system exerts the torque (note that u x , u y , and u z are the unit vectors of the body axes): 



The corresponding force transmitted at the center of mass of A is then: 


F* = 




and similarly for B: 


( 13 . 1 ) 
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F* + 


\ &B 


(13.2) 


In order to determine the constraint force F transmitted at the center of mass we write the 

equation of the variation of angular momentum for body A about its center of mass projected 
onto body axes: 


- [omM®/* + Va]»a - Mz X ¥ a + *A 


(14) 


where [I A ] is the inertia matrix of A in body axes andr^is the gravitational torque acting on A. 
After developing Eq. (14) and taking into account Eq. (1.1) we obtain: 


F = - 

r A.\ 


(I a x a I a . ' 

rr L+ -r- <w ,' 


and similarly for B: 

‘l n x B 1 


K — 


ft-**** 

b B °b 


Consequently, Eqs. (7) can be rewritten as follows: 

( ,.A 


*A + Met + [(w][a>]r A + 2[<o]r A - T*, z r A - 


\m A b A 2 m A b A ) m A 


tf±A_ , %d x A^ 

m A h A m A h A > 


U, 


(15.1) 


h + Me? + MMe? + 2[®]e? - r %.r B - 


l ,,B 


B ■ \ 


k e x B , U k B 

^ * 


\ m B^B 2 m BbB~ ) 


U, 


| h'B | h 

fibs" 


u v ~ 


el 

m B 


(15.2) 


m 4 m B 

r c —r 4 B -r B 


(15.3) 
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At this point we shall express Eqs. (15) in terms of the stretches x A , x B of the springs linking 
A with C and B with C. Let r A() ,r B0 ,r C0 be the radius vectors of the individual center of masses 
with respect to CM when the springs are at rest, then: 

r A = r A0 +L a u x = (r AOx +t A , r AOy , r AOz ) T (16.1) 

r B = r B0 +Cb U x = ( r BOx + Cb’ r B0y> f BOz) (16.2) 

r c = -m A /m c (r A0 +u A u x ) - m B /m c (r B0 +^ B u x ) (16.3) 

From Eq. (7.3) r c can be expressed as: 

r c =r co -l/m c ( m A ^ A + m B £ B ) u x (17) 

where 

r C o= -l/m c (m A r A() + m B r B0 ) (18) 

Let 8 A , 8 b be position vectors from the center of mass of C to the center of mass of A and B 
respectively, when the springs are at rest, then: 

8a= Tag- r co= r A o( 1 +m A /m c ) + r B0 m B /m c (19.1) 

&b= r Bo' r co= r Bo( 1 +m B /m c ) + r A0 m A /m c (19.2) 

We can now compute the stretches x A , x B of the two springs linking A with C and B with C 
as follows: 

X A = ( r A x -r A 0 x) - (r C x -r C ox) = Ca( 1 + m A /m c ) +C B m B /m c (20. 1 ) 

x B = (r Bx -r BOx ) - (r Cx -r a)x ) = ? B (1+ m B /m c ) +£ A m A /m c (20.2) 

Solving Eq. (16) for r A0 r B0 and Eq. (17) for £ A , £ B yields: 

r A o = 8A(m B +m c )/m T - 8 B m B /m T (21.1) 

r B0 = 8 B (m A +m c )/m T - 8 A m A /m T (2 1 .2) 

Ca = ^b niB^mr (21.3) 

C B = x B (m A +m c )/m T - x A m A /m T (2 1 .4) 

We can now write the equations of motions of the centers of mass in term of the variables 
x A (t) and x B (t) and the parameters 8 A , 8 B : 
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f„ + Wr, + [<»I«.]r. + W« - r|. : r A -( 

m A b A ) 
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i m A b A 2 m A b A 


m A 


* B + [»>„ + [*I4« ♦ w. - rU - ( %k 

\m B b B m B b B ) 
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( [b£b , h 

\m B b B 2 m B b B V j m B 
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m, m B 
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m r m r 


where 

r A = 8 A (m H +m c )/m T - 8 B m B /m x + (x A (m B +m c )/m T - x B m B /m T )u x 
r B = 8 B (m A +m c )/m T - 8 A m A /m T + (x B (m A +m c )/m x - x A m A /m T )u x 

After substituting Eqs. (23) into Eqs. (22), considering the x components of the Eqs. 
and (22.2) we obtain, for m c » m A and m B (as it is indeed the case) and - Fj* - 0: 


PA* A — (®v ■*" — "*■ ^ i x A* 2'*^ 

' ' m A b A~ m A b A 


4 S Az lj 7T I®' - 6 Av<»z ~ d Ax°>v ~ + &A <°x<*>z + t Av <O x U v 

\ m A b A/ • 


~~^Az.8xz ~ AySxy ^AxSxx ~ a x m ® 

PB%B -(«>v + «>? - 8xx) x B + —\ — x B + ~~~T X B 
' ' m B b B - m B b B 

+ S Bz ^-]®v - S B\<*>z - S Bx<Oy ~ &Bx<*>I + 6 B z<°x<»z + 6 Bv m x tO v 

\ m B b B ) ■ 

-&Bz8xz “ &By8xy + ^B.vSxv “ 0 
where 


( 22 . 1 ) 


( 22 . 2 ) 

(22.3) 

(23.1) 

(23.2) 

( 22 . 1 ) 


(24.1) 


(24.2) 
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m A^A '+- [a 


m B b B + I B 


Attitude equation of motion of the whole detector 

The attitude equations of motion can be derived based on the conservation of angular 
momentum. Let [I A ], [I B ], [l c ] be the inertia matrices of the three bodies in body axes <x,y,z. 
The total angular momentum of the system, projected upon <x,y,z> is written using the 
translation theorem of the angular momentum as follows: 

+ + [ / c] <B C +[fyfi]( Q Lt * + x + ®c x m c®c) ( 25 ) 

After some algebraic manipulations the derivative of the angular momentum projected onto 
body axes yield: 


— - Va]o>a + Vb]°>b + [/cK + [o>][ { a]o>a + [<»]Vb]o>b + H[/cK 

+m A r A x (f A + [d)]r A + [o)][<u]r A + 2[o>]f A ) 

+m B r B x (f fi + [d)]r s + [tu][(u]r B + 2[©]r B ) 

+m C r C * ( f C + W r C + MM r c + 2 [°>Yc) 

After taking into account Eqs. (22), Eq. (26) in body axes yields: 

[l A ](b A + [Ib]o>B + Vc]<»C + + [<*>][! b]°>B + 


m A + m A m c m A m B 

+ r A + r B 

1 m c m c 


< [rf v -]r A 


^0 X A _j_ Bd x _A u _ 

m A b A 2 m A b A 2 I X 


m kb A m A b A 


A ■ ] K 

, ®V K “ 


+ m B~ + m B m C + m A m B Ta 


lr k 0 X B , %0 X B „ h x B . l B 

1 a v: r B — t + — tt a rr + — T ~ 
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where x'™ is the gravitational torque on the whole instrument package. 

Numerical cases 


We integrated numerically the dynamics equations (both the attitude equations of the 
instrument package and the motion of the proof masses inside the detector) derived previously 
for roto-translational detector configuration. This configuration has appealing features from the 
point of view of constructing the suspensions of the proof masses and the ability of ground 
testing through gravitational excitation of the proof masses (see IFSI contribution Section). The 
new differential detector prototype was designed according to this configuration. 
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manufacture and provide an excellent frequency separation between their sensitive modal 
frequency (i.e., the one that bends the flexures) and any other modal frequency (see Ref. i). The 
acceleration detector prototype (see Figure 3) is shown here because it provides a good example 
of a roto-translational configuration. 

We explored cases with a precession frequency of the instrument package within the 0. 2-0.4 
Hz window that proved optimal for the translational detector (see Ref. i). As a general comment 
of simulation results, the roto-translational configuration has more internal dynamics than the 
translational configuration. The offsets between pivot axes and respective CMs of the proof 
masses couple the attitude dynamics of the instrument package to the internal dynamics of the 
proof masses more than with suspensions that are (ideally) located at the CMs of the proof 
masses. Moreover, the roto-translational configuration show the same sensitivity to gravity 
gradient forces (more on this point later) of the translational detector. This implies that the roto- 
translational configuration must satisfy the same tight requirement on the orientation of the spin 
axis of the translational configuration. After calling 0 the elevation angle of the spin axis with 
respect to the horizontal plane and 6 Z the centering error of the proof masses along the spin axis, 
then we must have 6 Z 0 < 1 micronxdeg. This requirement then drives the design of the leveling 
and release mechanism that needs to provide elevation angle errors less than 0.1 deg for 
centering errors of the proof masses equal to 10 microns. 

Unlike the translational detector configuration, the roto-translational configuration is also 
affected by inequalities between the moments of inertia about the pivot axes of the two proof 
masses. Unequal moments of inertia will produce slightly different rotations of the proof masses 
that are detected as a differential output at various frequencies 

A sample of simulations for a precession frequency of 0.2 Hz is shown in the following. 
Figure 4 presents the results for a baseline case with pivot offsets equal to 3.2 cm and relatively 
small centering errors of the proof masses. The output at low frequency is dominated by the 
precession dynamics-related noise that exhibits two peaks at the precession frequency and twice 
the precession frequency. Even after removal of those two harmonics by least-square fitting of 
the output in the time domain, the dynamics-related noise level at the signal frequency of 0.5 Hz 
is about 5x10 15 g. A better attenuation technique of the colored noise could improve the 
situation but, as shown later on, the dynamics-noise can be much more drastically reduced by 
changing the detector’s configuration. Note also the harmonic peak at 1 Hz which is associated 
with the gravity gradient forces due to the centering error of the proof masses. 

Subsequently, we have added an EP violation signal of 2x10 14 g to the previous case and the 
results are shown in Figure 5. The EP violation signal can be barely resolved after having 
filtered out the colored noise at low frequency as shown in the last panel of Figure 5. An EP 
violation signal of 5x10 14 g can be more clearly resolved but this is well above the sensitivity 
threshold level that we desire to achieve. The purely rotational configuration of the detector will 
cure this problem of high dynamics noise. 
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Damped 25% for 5 s, Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3, Iz = 0.1 8 (0.2 Hz prograde) 

OmegaX = 0.1 deg/s, Elevation angle = 0 deg; Spin = 0.5 Hz; Rotation Seq 2-1 -3, SmoothDamp = 0.04 
del-xA = del-yA= del-zA = 1 pm; del-xB = del-yB= del-zB = -1 pm; Pivots = 0.032cm 
Signal = 0 g; White Noise = 0 g; Initial xA = xB = +1 0 pm; Equations (8/1 7/04) 



Figure 4 Simulation results for a roto-translational configuration with instrument package’s 
precession frequency = 0.2 Hz and No EP violation signal. 
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Damped 25% for 5 s, Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3, Iz = 0.1 8 (0.2 Hz prograde) 

OmegaX = 0.1 deg/s, Elevation angle = 0.1 deg; Spin = 0.5 Hz; Rotation Seq 2-1-3, SmoothDamp = 0.04 
del-xA = del-yA= del-zA = 1 pm; del-xB = del-yB= del-zB = -1 pm; Pivots = 0.032cm 
Signal = 2x1 0 A -1 4 g; White Noise = 0 g; Initial xA = xB = ±1 0 pm; Equations (8/1 7/04) 



Figure 5 Simulation results for a roto-translational configuration with instrument package’s 
precession frequency = 0.2 Hz and an EP violation signal = 2x10 14 g. 
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Damped 25% for 5 s, Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3, Iz = 0.1 8 (0.2 Hz prograde) 

OmegaX = 0.1 deg/s, Elevation angle = 0.1 deg; Spin = 0.5 Hz; Rotation Seq 2-1-3, SmoothDamp = 0.04 
del-xA = del-yA= del-zA = 1 pm; del-xB = del-yB= del-zB = -1 pm; Pivots = 0.032cm 
Signal = 5x1 0 A -1 4 g; White Noise = 0 g; Initial xA = xB = ±1 0 pm; Equations (8/1 7/04) 



Figure 6 Simulation results for a roto-translational configuration with instrument package’s 
precession frequency = 0.2 Hz and an EP violation signal = 5x10 14 g. 
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Purely rotational configuration 

Detector configuration with pivot axes perpendicular to spin axis 

In a purely rotational configuration, the proof masses are constrained by the pivot axes 
(located at the CMs of the proof masses) to rotate with respect to the housing of the detector. A 
schematic of a configuration with the pivot axes of the proof masses perpendicular to the 
detector’s spin axis is shown in Figure 7. 


Pivot axis mass 2 



Figure 7 Schematic of rotational detector configurations with pivot axes of proof masses 
perpendicular to the detector’s spin axis: (a) 3-D view of one proof mass; and (b) geometrical 
arrangement of proof masses inside the detector. 

In this rotational configuration, one proof mass (labeled 1) is the “EPV-sensitive” mass that, 
thanks to its design consisting of two halves of different materials, is able to sense an EP 
violation force. Note also that an EPV force will produce a relatively strong torque (see later on) 
when compared to other perturbation torques thanks to its sizeable arm length with respect to the 
pivot axis. The second proof mass (labeled 2) is the “dynamical-reference” mass made of one of 
the materials of the first mass. As the name implies, the role of this proof mass is to move very 
closely (ideally in unison) to the other proof mass. This proof mass enables the removal of most 
of the instrument package dynamics from the detector’s output by differencing the individual 
motions of the two proof masses. 

Motion of proof masses with respect to detector’s housing 

Let us consider a set of three rigid bodies A, B, C constrained to one another in such a way 
that the relative motion of A with respect to C is a pure rotation around a pivot axis <dA> and the 
relative motion of B with respect to C is a pure rotation around a pivot axis <de> with <d A > and 
<d B > parallel to one another. 
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This time the pivot axes are ideally coincident with the CMs of the proof masses except for 
the construction errors. Let A and B be linked with C through a torsional spring dashpot system 
around the axes <d A > and <d B > characterized by elastic and damping coefficients (k A , § A ) and 
(k B , § B ) respectively. 



Figure 8 Relative position of individual centers of mass for pivot axis perpendicular to z-axis. 

Let <x,y,z> be a reference system centered at the global center of mass CM of the system 
composed by the three rigid bodies, which we assume coincident with the center of mass of body 
C, i.e., after having assumed that me » m A and m B . Let y be parallel to <d A > and <d B > (the y 
axis is exiting the page in Figure 8), z parallel to the central symmetry axis of body C, and x 
right-handed with y and z. A schematic of the relative positions of the centers of mass, inclusive 
of the construction errors 6 A and 6 B , is shown in Figure 8. 

The tangential unit vectors Uj A and Utb can be written as: 

u ta ■ j— ,-6 az) T 

°Axz 

u tb ■ t — ( 6 fiv ,0,-6 fi -) r 

°Bx- 

with 

&Axz “ Alv 2 + &A : 2 '■> &Bxz = 
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l 


Let cdc = (co x ,co y ,a) z ) T be the angular velocity vector of body C and, consequently, of the 
reference frame <x,y,z>. The angular velocities of body A and B, projected on <x,y,z> are: 


®A 


\T 

+ -T 4 -® J 

°Axz ) 


a> B 




Vb 


6 


,© 7 


Bxz 


(28.1) 


(28.2) 


where r\ A and q B are the tangential displacements of the center of masses of A and B. with 
respect to the initial geometry conditions drawn in Fig. 8. The equations of motion of the system 
written with respect to an inertial system centered at the Earth’s center are: 


m A ^CM + p A )' - FJ + F k A + F? + fJ + m A • af ’ v 

(29.1) 

m B^CM +Pfl)" F J + + f| + Fg 

(29.2) 


where p is the radius vector of each individual center of mass relative to CM and projected on 
the inertial system, R CM is the radius vector of CM from the earth center projected on the inertial 
system, F g is the earth gravity field force, F k the equivalent force of the internal spring and F 5 the 
equivalent force of the internal damper acting on each mass. F b is the force transmitted to each 
mass from the pivotal constraint as a consequence of the rigid body attitude motion of A and B. 
We also added the acceleration a EPV acting on the mass A to simulate a violation of the 
equivalent principle. The equations of motion of the system center of mass are: 


R 


CM 


H e R 


E IV CM 
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3 

CM 



(30) 


where F c . is a perturbing force acting on the external case during free fall and m x the total mass 
of the system. The gravitational forces, appearing in Eqs. (29), can be expressed like Eqs. (4) of 
the previous subsection. 


After plugging Eqs. (30) into Eqs. (29) we obtain: 
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(31.1) 
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The previous equations can be written with respect to body axes <x,y,z> as follows: 
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(32.1) 

(32.2) 


where r A ,r B are the radius vectors of each individual center of mass relative to CM and projected 
on the body frame; R IB , R BI are the rotation matrixes from the inertial to the body frame and 
viceversa, [w] is angular velocity matrix of the body frame, and is the Earth’s 

gravity gradient matrix projected in body frame (see previous subsection for details). The forces 
F* and F' transmitted by the rotational spring dash-pot system can be obtained by following the 
same procedure outlined in the previous subsection. 

A rotation of the proof mass A, «■ Va/^a.xz w * t * 1 respect to C forces the spring dash-pot 
system exerts the torque: 

t sd ■ (-*0 1 &a -£e u y (33) 


The corresponding force transmitted at the center of mass of A is then: 



&Axz &Axz t 


«7A 


(34) 


and similarly for B with 0 B - tj b /8 Bxz : 
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In order to determine the constraint force F ft transmitted at the center of mass we write down 

the equation of the variation of angular momentum for body A about its center of mass projected 
onto body axes: 

" MM® a + M® a - I-A * Fa + *A (35) 


where [I A ] is the inertia matrix of A in body axes and x g A is the gravitational torque acting on A. 
Consequently, after developing Eq. (35) and taking into account Eq. (28.1) we obtain: 
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and similarly for B: 


t B 


V d Bx: °Bxz ) 


Hence Eqs. (32) can be rewritten as follows: 
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(36.2) 


(36.3) 


At this point we can express Eqs. (36) in terms of the stretches x A , x B of the springs linking A 
and B with C, as follows: 


r A = &a +*1 a Uta =(S Ax + 0 a 6 Az ,0,6az ~ 0 a&Ax) T 

r B = 6 b +t 1bU T b = (&BX + 0 B&B-’°' 6 Bz “ 0 B 0 Bxf 


(37.1) 

(37.2) 


Finally, projecting Eq. (36.1) of proof mass A along the tangential direction u TA and dividing 
by the arm length 5 Axz in order to express the equation in terms of the relative rotation 0 A : 


0 a + C 2 6 a +C>6 a + (b 


+C 4 (cb x + (O y O) z - g v .) + C 6 (w. + (o x io y + g, v ) 


(38) 
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(39) 


with 1 Ay m I Ay + m A&Axz *^Arz " 'J^Ax ^Az • 

The equation relative to the rotation 0 B of mass B are written by simply replacing the 
subscript A with B and setting the equivalence violation torque to zero. Note that the assumption 
of m c » m A ,m B eliminates the cross-coupling from the two equations of motion. 

Effect of gravity gradient torques 

The attitude equations of motion are not shown for the purely rotational case because their 
derivation follows the same strategy outlined previously for the roto-translational detector. 
However, there are important points to be made about the effects of the gravity gradient torque 
acting on the instrument package versus the Equivalence Principle violation (EPV) torque. The 
components in body axes of the gravity gradient torque, for the 2-1-3 Euler’s rotation sequence, 
are: 


V GG - -3 — - / v )cosV>cos0sin0 


X GG - -3^p(/ x - 7.)sin^cos0sin0 


R 

gg o GM 

If -7 ■ 


3 — — (/ v - I^sinrjfcosrjfcos 2 0 


(40) 


where l x = / A( + I Bx + I Cx , I y = I Ay + I By + I Cy and I. = I Az + I B . + l Cz are the overall moments of 
inertia of the instrument package. 

The component of interest of the EPV torque is the one that produces a differential motion 
between the proof masses, i.e., the component about the pivot axis that we call generically as the 
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v-axis with v = x, y, or z. The EPV torque is far too small to produce any effects about the other 
axes. After calling p the arm length of the EPV violation acceleration vector with respect to the 
pivot axis, the v-axis component of the EPV torque can be computed as: 


EPV 

T v 



(41) 


EPV • 

where the EPV acceleration in body axes is obtained by transforming the inertial-axes 
EPV 

acceleration a as follows: 

„EPV r j? 1» EPV 

The expression of the sensitive EPV torque component (i.e., about the pivot axis) changes 
only slightly depending on the selection of the pivot axis-v alignment with any of the body axes, 
as follows: 


x EPV «_p_ cos ^; cos 0 \j EPV 

Xy PV m-p z smiltcosOa EPV (42) 

x EPV m p y s\MltcosQa EPV 

In Eqs. (42), the spin angle (about the z-axis) rp = ul (where co is the spin rate), 0 is the 
elevation angle, and <j)the azimuth angle. The key observation here is that the EP violation 
torque is always modulated at the frequency ou, irrespective of the orientation of the pivot axis 
along any one of the body axes. 

On the contrary, the body-axes components of the gravity gradient torque (of Eqs. 40) have 
different frequency content. Specifically, the components about the x and y axes are modulated 
at ou while the component along the z-axis (i.e., the spin axis) is modulated at 2to. This 
consideration has very important consequences for acceleration detectors that utilize rotational 
configurations. The issue here is the gravity gradient torque acting on the instrument package 
(see Eqs. 40) because the instrument package needs to have for dynamic reasons a longitudinal 
moment of inertia different from the transverse moment of inertia and, consequently, this gravity 
gradient torque is significant. Gravity gradient torques acting on the proof masses are much 
smaller than on the instrument package because the inertia ellipsoids of the proof masses are as 
spherical as possible (within the construction tolerances). 

Detector configuration with pivot axes parallel to spin axis 

In order to separate the effect of the gravity gradient torque acting on the instrument package 
from an EP violation signal, it is convenient to align the pivot axes of the proof masses with the 
symmetry axis z of the instrument package, as shown in Fig. 9. 
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Figure 9 Schematic of purely-rotational detector configuration with pivot axes of proof masses 
parallel to the detector’s spin axis: (a) 3-D sketch of proof mass; and (b) geometrical 
arrangement of proof masses inside the detector. 

The equations of motion of the proof masses for this configuration can be derived from Eqs. 
(38) and (39) through simple axes transformations to obtain: 

d A + E 2 6 a + E y d A + 6). 
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Note that for the purely rotational configuration the EPV torque per unit mass is as follows: 

= p x a £P ' / cos(tur)cos(0) (45) 


where p x , the arm length, is of order a few centimeters because half of one proof mass is made of 
a material different from the other half. The size of the arm length has very positive 
consequences regarding the effects of the Earth’s gravity gradient and ultimately the tolerance to 
attitude errors of the spin axis. 

An inspection of Eqs. (43) and (44) shows that the gravity gradient “force-related” torques 

GG 

per unit mass (i.e., excluding the “moment-related” gravity gradient torque t A; ) are all equal to: 
l 1- !: I = ^Ai^AjShk with i,j,h,k = X,y,z. (46) 


The “force-related” gravity gradient terms of greatest concern are the off-diagonal ones (with 
h * k) that are modulated at the spin frequency w. The most important among those are the g yz 
and g xz terms that have the form: 


8zx 


8xz 



sin(cu/) cos(0)sin(0) 


8 Z y 


gyz 



cos (iot) cos(0)sin(0) 


These terms, which are directly proportional (for small angles) to the elevation angle 0 (while 
the g xy term is not), drive the requirement on the tolerable attitude error of the spin axis with 
respect to the horizontal plane. We then compare the maximum value of the EPV torque of Eq. 
(45) to that of the gravity gradient torque of Eq. (46) for typical values of the differential 
acceleration detector as follows: a EPV = 5xl0 14 m/s 2 , p x = 3 cm, R = 6418 km. Finally, we obtain 
that for very conservative values of differential errors (i.e., between the two proof masses) 6 X , 6 y 
and b 7 < 50 microns and an attitude error of the spin axis < 7 deg, the (co-modulated) gravity 
gradient torques are smaller than the EPV torque at the threshold sensitivity of the detector. In 
other words, this design of the differential acceleration detector has cured the problem of the 
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tight requirement on the orientation of the spin axis during free fall and also the centering errors 
between the proof masses CMs. As a result, the requirement on the horizontality of the spin axis 
at release can be relaxed, thus simplifying the design of the leveling mechanism and the 
spin/release mechanism. 


The orientation of the pivot axis along the symmetry axis z also cures the problems generated 
by the “moment-related” gravity gradient torques r\ z and r|, acting on the proof masses 
because the z-axis components are now modulated at 2a> [see the third of Eqs. (40)] and, 
consequently, their frequency is well separated from the EPV signal. This is also good news 
because a very strict equality between the moments of inertia of the two proof masses is no 
longer required. However, the ellipsoid of inertia of each proof mass needs to be spherical with 
6L/L < 10 4 (in which L is a linear dimension of the proof mass) in order to make negligible the 
higher-order mass moments which are coupled to nearby masses as explained in Ref. i. This 
issue is common to all configurations of the detector. 


In other words , the purely rotational configuration with pivot axes aligned along the spin 
axis of the detector has mitigated, thanks to frequency separation, the issues associated with 


gravity gradient torques acting on both the proof masses and the instrument package. Moreover, 
the combination of the torsional degrees of freedom and the design of the ”EPV-sensitive” proof 
mass made of two different materials has maximized the effect of an EP violation signal, that 
now overpowers the effects of the gravity gradient forces. 


Numerical cases 

We ran several simulations of the two configurations of the purely rotational detector: (a) 
with the pivot axis orthogonal to the spin axis; and (b) the pivot axis parallel to the spin axis. 
Note that in all those simulations the proof masses are initialized with very large oscillation 
amplitudes (up to their full range), which we expect to occur after detector’s release. Those 
oscillations are abated during the first few seconds by the feedback loops (see Ref. i), the 
dynamics effects of which are incorporated into the simulations. At the 5-s mark, the feedback 
loops are turned off and the detector enter its high-Q (i.e., negligible damping) mode. 

The detector is also released with transverse (orthogonal to the spin axis) rotational velocity 
errors of 0. 1 deg/s that is the maximum error expected from the spin/release mechanism. An 
inclination error of the spin axis at release is also considered, ranging from 0. 1 deg for the earlier 
simulations to 1 deg for other simulations. 

Pivot axis perpendicular to spin axis 

First, we adopt a precession frequency for the instrument package of 0.3 Hz which is in the 
previously-defined desirable frequency window. Besides the typical errors at release, we also 
assume an error between the elastic oscillation frequencies of the proof masses 6to/u) = I0 \ 
which is consistent with a common-mode rejection factor (CMRF) better than 10" 4 . The CMRF 
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depends on the equality of the elastic frequencies of the proof masses. These two frequencies 
can be independently calibrated by adjusting the feedback loops of the proof masses. 

Figure 10 shows the results of a simulation without any EP violation signal. Note that there 
are no significant components appearing at the signal frequency of 0.5 Hz (i.e., at the spin 
frequency) and that there is also no significant harmonic component at twice the precession 
frequency. Figure 1 1 shows simulation results of a similar case in which we have added an EP 
violation signal of 1x10 14 g. The original signal can be recovered with enough precision by 
means of frequency analysis. 

We need to remind that for proof masses with centering errors, the precession dynamics of 
the package generates differential acceleration harmonics at the spin frequency and twice the 
spin frequency (see Ref. i). A notable advantage of the purely rotational configuration is that the 
amplitudes of both harmonics are attenuated and the harmonic at twice the spin frequency 
becomes negligible with respect to other noise sources. Consequently, we are freer than in the 
other configurations in choosing the value of precession frequency because the harmonics higher 
than the precession frequency will not significantly affect the signal frequency. The final result 
is that it is convenient to lower the package’s precession frequency in order to reduce the 
amplitude of the precession-dynamics-associated noise. Since we would like to resolve the 
precession dynamics out of the output data through frequency analysis, it makes sense to choose 
a precession frequency such that at least a couple of full precession cycles occur during free fall. 
In conclusion, a precession frequency of 0.1 Hz appears an optimal value for a free fall duration 
of about 25 s. 

Figure 12 illustrates the point that, if the moments of inertia of the two proof masses are 
perfectly identical, then the noise level at the signal frequency is low. An EPV violation signal 
(here assumed equal to 1x10 14 g) would be clearly detectable above noise as shown in Fig. 12. 
However, if we now assume a fractional error of 2x10 3 (i.e., an absolute error of 10 5 kg-m 2 ) 
between the moments of inertia of the proof masses, the gravity gradient torque associated with 
that error, which appears at the same frequency of the signal, overpowers the signal as shown in 
Fig. 13. 

In other words, the purely rotational configuration with the pivot axis perpendicular to the 
spin axis has mitigated the impact of the precession dynamics on the accelerometer output but 
has not improved the issue of “moment-related” gravity gradient torques acting on the proof 
masses themselves and the effect of an inequality in moments of inertia on the differential 
acceleration. 
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TORSION DETECTOR; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.1 2 (0.3 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 0.1 deg; Spin = 0.5 Hz; lx * ly; OM_err = 10 A -5 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; Damp = 37% 

Signal = 0 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations (1 2/05/04) 



Figure 10 Simulation results for the rotational configuration with pivot axes orthogonal to the 
spin axis, instrument package precession frequency = 0.3 Hz, a difference between oscillation 
frequencies of the proof masses 6a>/oa = 10' 5 , and no EP violation signal. 
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ROTATION DETECTOR; Elastic freq. = 3 Hz; lx = 0.3003; ly = 0.3, Iz = 0.12 (0.3 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_Err = 10 A -5 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; Smooth damping = 25% 
Signal = 1x1 0 A -1 4 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations (12/07/04) 



Figure 11 Simulation results for the rotational configuration with pivot axes orthogonal to the 
spin axis , instrument package precession frequency = 0.3 Hz, a difference between oscillation 
frequencies of the proof masses 6u>/u> = 10'\ and an EP violation signal = 1x10 14 g. 
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TORSION DETECTOR; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 0.1 deg; Spin = 0.5 Hz; lx * ly; lAy = IBy 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; 5-s, Damp = 37% 

Signal = 0 g; Noise = 0 g; Initial xA = xB = ±60 pm; Equations (12/05/04) 



Figure 12 Simulation results for the rotational configuration with pivot axes orthogonal to the 
spin axis, instrument package’s precession frequency = 0.1 Hz, perfectly spherical ellipsoids of 
inertia of the proof masses, and no EP violation signal. 
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TORSION DETECTOR; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 0.1 deg; Spin = 0.5 Hz; lx * ly; lAy = IBy 
del-xA = del-yA = del-zA = 1pm; del-xB = del-yB = del-zB = -1pm; 5-s, damp = 37% 

Signal = 1 xl 0 A -1 4 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations (12/05/04) 



Figure 13 Simulation results for the rotational configuration with pivot axes orthogonal to the 
spin axis, instrument package precession frequency = 0.1 Hz, perfectly spherical ellipsoids of 
inertia of the proof masses, and an EP violation signal = 1x10 14 g, which is accurately recovered 
by the frequency analysis. 
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TORSION DETECTOR; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle - 0.1 deg; Spin = 0.5 Hz; lx * ly; OM_FErr = 5x10 A -6 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; IBy&IBz_FErr= 2x1 0 A -3 
Signal = 1x1 0 A -1 4 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations ( 1 2/05/04) 



Figure 14 Simulation results for the rotational configuration with pivot axes orthogonal to the 
spin axis , instrument package’s precession frequency = 0.1 Hz, EP violation signal = 1x10 14 g 
and a fractional error between the moment of inertia of the proof masses of 2x1 0' 3 . Note that the 
error associated with the moments of inertia has overpowered the EP violation signal. 
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Pivot axis parallel to spin axis 

As mentioned previously, the configuration with the pivot axes of the proof masses parallel 
to the spin axis of the detector cures the issue of the gravity gradient torques. In this geometrical 
arrangement, the effects of those torques appear in the accelerometer output at twice the signal 
frequency and, consequently, do not overlap with the signal. 

In the first simulation (shown in Fig. 15), we have assumed the usual velocity error at release 
of 0. 1 deg/s (about an axis perpendicular to the spin axis) and we have increased the error on the 
elevation angle of the spin axis to 1 deg. The precession frequency of the package is 0. 1 Hz. 
The fractional error between the proof masses elastic frequencies is 6co/o) = 5x10 6 (which 
provides a common-mode rejection factor CMRF » 10 5 ), the fractional error between the proof 
masses moments of inertia 51/1 = 10 5 , and no EP violation signal. The figure clearly shows that 
the level of dynamics-related noise at the signal frequency of 0.5 Hz is very low. 

Figure 16 shows a case with same error parameters as in the previous simulation but with the 
addition of an EP violation signal = 1x10 14 g. In the simulation of Fig. 17, we have increased 
the fractional error between moments of inertia of the proof masses to a very lax 10 3 in order to 
show that, thanks to frequency separation, the gravity gradient torques (of the second-order) are 
not an issue with this detector’s configuration. Actually, the tolerable difference in inertia 
characteristics between the proof masses will be limited to a stricter value by the higher-order 
moments coupled to nearby masses. The higher-order mass moments require shape accuracies 
between the two masses better than 6L/L = 10 4 which in turns will imply second-order moments 
of inertia with differential errors 6I/I < 10^. 

In the simulation run of Fig. 18, we have increased the CM centering errors further and 
relaxed the difference in natural frequencies between proof masses as follows: 6 A = 25 microns, 
6 b = 30 microns 6u)/u) = 10' 5 , 6I/I = 10 4 . Moreover, we have reduced the EP violation signal to 
2x10 15 g. As shown in Fig. 18, the small signal has been successfully resolved from the 
dynamics-related noise and gravity gradients. In other words, the detector sensitivity is no 
longer limited by the dynamics-related noise and gravity gradients but rather by the intrinsic 
noise of the differential accelerometer. The latter noise can be modeled as white Gaussian noise, 
which was assumed null in these simulations (see later on for the extraction of signal from noise) 
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TORSION DETECTOR; PIVOT-z; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_FErr = 5x10 A -6 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; IBx_FErr = 1 0 A -5 
Signal = 0 g; White Noise = 0 g; Initial xA = xB = +50 pm; Equations (01/03/05) 



Figure 15 Simulation results for the rotational configuration with pivot axes parallel to spin 
axis, instrument package’s precession frequency = 0.1 Hz, proof masses inertia error 51/1= 10 5 , 
oscillation frequency error 6u)/u> = 5xl0 6 , and no EP violation signal. 
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TORSION DETECTOR; PIVOT-z; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_FErr = 5x10 A -6 
del-xA = del-yA = del-zA = 1 Oym; del-xB = del-yB = del-zB = 1 Zfjm; IBx_FErr = 1 0 A -5 
Signal = lx!0 A -14 g; White Noise - 0 g; Initial xA = xB = ±50 (jm ; Equations (01/03/05) 



Figure 16 Simulation results for the rotational configuration with pivot axes parallel to spin 
axis , same errors as in Fig. 15, and an EP violation signal = 1x10 14 g. 
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TORSION DETECTOR; PIVOT-z; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_FErr = 5xlO A -6 
del-xA = del-yA = del-zA = 1 0pm; del-xB = del-yB = del-zB = 1 2pm; IBx_FErr = 1 0 A -5 
Signal =1x1 0 A -1 4 g; White Noise = 0 g; Initial xA = xB = ±50 pm; Equations (01/03/05) 



Figure 17 Simulation results for the rotational configuration with pivot axes parallel to spin 
axis , instrument package’s precession frequency = 0.1 Hz, proof masses inertia error 51/1 = 10 \ 
oscillation frequency error 6(o/co = 5x1 0'\ EP violation signal = lxlO 14 g. 
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TORSION DETECTOR; PIVOT-z; Elastic freq. = 3 Hz; lx - 0.3; ly = 0.3003, Iz = 0.24 (0.1 Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_FErr - 5x10 A -6 
del-xA = del-yA = del-zA = 1 Opm; del-xB = del-yB = del-zB = 1 2pm; IBx_FErr = 1 0 A -3 
Signal = 1 xl 0 A -1 4 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations (01 /03/05) 



Figure 18 Simulation results for the rotational configuration with pivot axes parallel to spin 
axis, instrument package’s precession frequency = 0.1 Hz, proof masses inertia error dl/I = 10'\ 
oscillation frequency error 6w/u) = 5x10 \ EP violation signal = 1x10 14 g. 
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TORSION DETECTOR; PIVOT-z; Elastic freq. = 3 Hz; lx = 0.3; ly = 0.3003, iz = 0.24 (0.1Hz prograde) 
OmegaX = 0.1 deg/s, Elevation angle = 1 deg; Spin = 0.5 Hz; OM_FErr = 1 0 A -5 
del-xA * del-yA = del-zA = 25pm; del-xB = del-yB = del-zB = 30pm; IBz _FErr = 1 0 A -4 
Signal = 2x1 0 A -1 5 g; Noise = 0 g; Initial xA = xB = ±50 pm; Equations (01 /03/05) 



Figure 19 Simulation results for rotational configuration with pivot axes parallel to spin axis, 
instrument package’s precession frequency = 0.1 Hz, proof masses inertia error bill = 10' 4 , 
oscillation frequency error Sco/co = 10‘ 5 , large centering errors of CMs, spin axis elevation error = 
1°, and EP violation signal = 2x10 15 g. No white noise was added to the accelerometer output. 
Note that the small signal is successfully resolved out of the dynamics-related noise. 
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Detector conceptual design 


Proof mass geometry 

The analysis conducted thus far has resulted in a conceptual design of the flight-experiment 
detector that is able to reject almost any reasonable perturbation and, at the same time, magnify 
the effect of an EP violation signal with respect to perturbations. A large number of numerical 
simulations carried out for the various configurations have proved that the purely-rotational 
configuration with pivot axes parallel to the detector’s spin axis can effectively separate the 
dynamics-related effects and gravity gradients from the signal. 


An important element of the rotational detector is the design of the proof masses, consisting 
of: (1) an EPV-sensitive mass made of two different materials; and (2) a dynamics-reference 
mass that enables the attenuation of the rigid-body dynamics and its effect on the detector’s 
differential output. In this section we would like to dwell further on the inertia characteristics of 
the proof masses and options for their mechanical design. 



Figure 20 Geometry of proof masses. 

Equal masses of the proof masses are not required because, thanks to the EP, the mass does 
not count in free fall. The equality of the oscillation frequencies of the proof masses is important 
in a differential detector because we want the two masses move in synchronicity and with equal 
amplitude under the effect of a common-mode excitation. The frequencies depend on the ratio 
K/I where K is the rotational elastic stiffness and / the moment of inertia of each mass about the 
pivot axes. Although in general we could design for equal ratios, here we focus on the case in 
which the stiffnesses are the same and, consequently, we need to design for I A = /„. Remember 
also that each proof mass needs to have a spherical ellipsoid of inertia which drives the design 
into using cylinders with a specific ratio between radius and length as shown later on. 
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In the schematic drawing of Fig. 19 we use the subscript 1 for the internal and 2 for the 
external proof mass; l, and l 2 are the cylinders’ lengths, R el and R e2 are the external radii and R, , 
and R i2 the internal radii of the proof masses. Moreover, we adopt a body-attached reference 
frame with z along the longitudinal axis of the proof mass and x and y along the transverse axes. 
Note that the naming of the axes is conventional and, consequently, the analysis shown here is 
valid for either a configuration with the pivot axis perpendicular to the spin axis or the pivot axis 
parallel to the spin axis. 


We need to design for proof masses with: (a) spherical inertia ellipsoids (i.e., same moments 
of inertia about the three axes); (b) same moments of inertia between one another; (c) and 
coincident centers of mass in the geometrical centers, where the pivots are attached. Let us focus 
first on the internal proof mass which consists of a cylinder with two halves (named a and b) 
made of different materials. We first define the moments of inertia about the three axes as 
follows: 


rCM 

hz 




,CM 

1 


'1 ax 


+ I lbx + (™l« + '"!*)( 


jCM ,CM 
Mv "'I* 


(47) 


where the asterisk indicates that that moment of inertia is computed with respect to the center of 
mass of the half portion of the proof mass while the superscript CM stands for the overall center 
of mass. 


The condition of spherical inertia ellipsoid that is, I lx = /,. = 7, v , can be expressed as: 


'laz 


hbz " hax + I \bx+ ( m l a + 


(48) 


in which the cylindrical symmetry makes the second equality redundant. The condition of CM in 
the middle of the proof mass implies obviously that 


1 

m \a mfn \b m 2 m l 


(49) 


After expressing the moments of inertia in terms of geometrical parameters, we obtain 
haz m ^ m \[ R e\a + R f\aj (50.1) 

hbz-^ m \{ R e\b* R iu) <50.2) 
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hax ™ 24 - > Ri\ a + (/] 12) J 
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(50.3) 
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hbx-^rn, 

Since we have 5 geometrical variables and one equation, we need to impose more conditions 

in order to restrict the numbers of possible solutions. If we require that the two halves of the 

* * 

internal proof mass have equal moments of inertia about their own CMs then l\ az — will 
yield 


(4« + RL)/2-(A 2 i(. + R®)/2 


(51) 


Equation (51) restricts the possible solutions to only those pairs of radii with equal 
geometrical means for the two halves of the proof mass. 


If we now impose that - I^ x then Eqs. (47) and (48) yield: 
2 /n i(^ 2 i +R n) 
f x M 3R?i + if] 


,CM 1 
7,5 2 


(52.1) 

(52.2) 


which, when equated to one another, yields the usual relationship for cylinders with spherical 


inertia ellipsoids 



(53) 

After calling p a and p b the densities of the materials that make up the two halves (in which 
we have assumed p a < p b ), we have 

m \ a m Pa^^ R e\a ~ R i 1 a J 

(54.1) 

m \b " PbX^elb ~ 

(54.2) 

where m, = m la + m ]b . Finally, Eqs. (49), (51) and (54) lead to: 


R el a " + Pa! Pb ) R i 1 /) 0 “ Pa/ 1 Pb) j 

(55.1) 
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R ha m Pa / Pb) + R hb{ [ + Pa /Pfc)] 


(55.2) 


which relate the radii of one side to the other side of the proof mass. 

The relevant equations for the external proof mass (subscript = 2) can be simply summarized 
as follows: 


m 2-Pa*( R e2- R h)l2 

jCM _ jCM jCM 1 ( p2 . n2 \ 

l 2x l 2y ml 2z " 2 m 2 \ R e 2 + R i2) 

h “ V3V4 + 


(56.1) 

(56.2) 

(56.3) 


where in Eq. (56.1) we have assumed that the density of proof mass 2 is equal to the density of 
one of the halves of proof mass 1. The previous equations can be used to dimension the two 
proof masses so that they satisfy the conditions imposed on their inertial characteristics. 



Figure 21 Geometry of internal proof mass with flat inner and outer surfaces. 


In a further refinement of the internal proof mass design, we also impose that the inner and 
outer radii are the same for the two halves. These additional geometrical constraints will 
simplify the design of the signal pick up plates of the accelerometer because, in this case, the 
internal proof mass will have flat inner and outer surfaces. However, in order to satisfy these 
two additional constraints without violating the previous ones, we have to split the right half of 
the proof mass into two concentric cylinders: the inner one will be mounted flush with the inner 
radius of the left half and, likewise, the outer one will be mounted flush with the outer radius. In 
other words, the inner and outer radii of the left half have been replaced by the inner radius of the 
outer cylinder of the outer radius of the inner cylinder (see Fig. 21). By following the new 

parameterization, the conditions l\ ax - l\ bx and l\ a: - l\ bz , combined with Eqs. (54), yield: 
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-(l-Pa/PtK 



(57.1) 

(57.2) 


in which we have dropped the subscript 1 for the sake of simplicity. Equations (57) replaces 
Eqs. (55) for the flush-cylinder configuration. 


In a modification of the previous design, the internal proof mass would involve the inner 
cylinder of the right half to consist of two layers: the innermost one (on the side of the pivot 
attachments) made of the same material of the left half and followed by a thicker cylinder of the 
higher-density material (see Fig. 22). The outer cylinder will also be made of the higher-density 
material. The goal obviously is to maximize the amount of mass made of the higher-density 
material while at the same time provide for ease of manufacturing and construction precision. 
Having the innermost surface of the internal proof mass made of the same material could help in 
constructing the pivots of the proof masses which are also made of the same materiaf\ Other 
geometrical options are also possible. We have not explored them further because at this point 
of the project we are more interested in developing the conceptual models and analyzing their 
merits rather than focusing on the detailed engineering of the flight detector. The latter effort 
will require a thorough design and development work of the flight detector to be led by our 
partners at the IFSI laboratory. 


Pivots 



Figure 22 Schematic of internal proof mass with inside surface of cylinder made of one material. 
Sample dimensioning of proof masses 

We have utilized the previous equations to derive a sample design of the two proof masses tat 
could be used in a purely rotational configuration of the detector. This is only a sample design of 
the proof masses that does not intend to freeze the dimensioning of the actual flight detector for 
the EP experiments. In the example of Figure 23, we have adopted Aluminum as the low-density 
material and a Copper alloy as the high-density material. Aluminum and Copper are materials 
that have been widely used for EP testing in the past and they are adopted here for 
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exemplification only. Other materials pairs, like Beryllium and Platinum-Iridium for example, 
would lead to a design of the proof masses with different dimensions. 



«*— * 


Figure 23 Sample dimensioning of proof masses made of Aluminum and Copper-alloy. 

Many choices of materials are possible because the equations derived above have solutions 
for any ratio of the material densities p^p,,. In the design shown in Fig. 23, the inertia 
characteristics of the proof masses are as follows: 

m, = 2.82 kg; I xl = I yl = I zl = 5.47x10 3 kg-m 2 (1 = internal proof mass) 

m, = 1.12 kg; I x2 = I y2 = I z2 = 5.47x10 3 kg-m 2 (2 = external proof mass) 

These values have been used in the dynamics simulations shown in the previous section for 
the rotational configuration of the detector. As stated previously, this dimensioning is one 
possibility among many options. 

Concluding remarks 

The analysis of the dynamics in free fall of various detector configurations in the presence of 
release errors and proof masses imperfections has led to the design of a suitable detector 
configuration. This configuration provides a very low level of dynamics-related noise at the 
signal frequency that meets the desired threshold sensitivity goal for the detector. The 
configuration consists of proof masses that have purely rotational degrees of freedom. The 
internal proof mass is the “EP-violation-sensitive” mass that, thanks to its design made of two 
halves of different materials, is able to sense an EP violation force. The second mass, made of 
one of the materials of the other mass, is the “dynamical-reference” mass. The role of this proof 
mass is to attenuate (and ideally remove), through output differencing, the dynamics-related 
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noise from the detector’s output. Moreover, with this design an EP violation force will produce a 
relatively strong torque when compared to other perturbation torques. 

We have demonstrated through extensive simulations that this detector design is capable of 
detecting a very small differential signal and resolving that signal out of the dynamics-related 
noise and gravity gradients in the presence of conservative values of the construction errors of 
the proof masses. This conclusion implies that, with this new design, the dynamics-related noise 
is no longer the limiting factor of the detector’s threshold sensitivity’. The limiting factor becomes 
the intrinsic noise of the detector (e.g., the Brownian noise of the proof masses). 

We have also derived the geometrical relationships between the internal and external proof 
masses that provide the desired inertia ellipsoids for each mass and the equality between the 
moments of inertia of the two proof masses. Moreover, we have worked out design options for 
the internal proof mass (the one made of two different materials) and the geometries that satisfy 
the inertia and mass requirements for that proof mass. 
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Detection and Estimation of the EP Violation Signal 

Characteristics of the data 

The sampled data consists of a colored noise, a white Gaussian noise (WGN) and a possible 
EP violation signal. The colored noise is produced by the precession dynamics (at a known 
frequency) and by elastic modes of vibration. The Signal to Noise Ratio due to the colored noise 
is very small and could be as much as -80dB. The white noise represents a thermal noise and 
may be at the same order of magnitude as the violation signal. The violation signal is at a known 
frequency (0.5 Hz) and phase. In this report we will call the EP violation signal the signal. 
Accordingly, there are two major difficulties in extracting the signal. The first problem is the low 
SNR and the possible masking of the signal by colored noise components (e.g., the harmonics 
associated with the precession dynamics). The second problem is the extraction of the signal 
from the white noise. Since the harmonic related to the signal frequency may be affected by the 
white noise, we need to estimate the confidence of detection as well. For that purpose we will 
assume that there is no component of colored noise at the signal frequency. 

The data record length is T - 20 seconds and the sampling rate is typically f s - 100 Hz. 
Therefore the data consists of 2000 samples. Accordingly, the spectral resolution is 
A/'- l/Tm 0.05 Hz. The bin width on the other hand is Vw, - f s /N dft where N DFT is the 
number of points in the DFT. The maximum number of point effectively used for an FFT (out of 
a sample of 2000 ) is 1024 that gives a bin width of about 0.097 Hz. Therefore we prefer a DFT 
algorithms of 2000 points, with a bin width of 0.05 Hz, rather than an FFT algorithms. In 
addition to sufficient resolution we require a bin at 0.5 Hz (the assumed signal frequency, /q). 
The condition for that is the existence of an integer K such that K - N f$l f s . Therefore 
N ~200K for the current data. The conclusion from this discussion is that an acceptable 
N dft for bin widths less than 0.1 Hz are 1000, 1200, 1400, 1600, 1800 and 2000 for a bin width 
of 0.05 Hz. 

As a test-bed we used four sets of sampled data associated with different cases and/or 
detector configurations. The data include a colored noise and an assumed signal of 2.8xl0 16 m 
(equivalent to a differential acceleration of 1x10 14 g). The sampled data and their DFTs are 
shown in Figs. 24 and 25. It is obvious that the DFT alone is not sufficient to resolve the signal. 
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Figure 25 DFT of sampled data without pre-processing. 

Our goal was to search and test feasible methods for detecting and estimating the signal. In 
this report we deal separately with the colored noise and with the white noise. We first present 
few methods to pre-process the data in order to improve the signal detection. We then test few 
filtering techniques to estimate the signal with the pre-processed colored noise. Lastly, we 
assume a pure white noise and we estimate and detect the signal by using an optimal statistical 
filter. 

Data pre-processing 
Windowing 

In order to enhance the observability (frequency detection) of the DFT data we need to pre- 
process the data through a window. We tested three window algorithms: 

Hamming: W (n) ■ 0.54 - 0.46cos(2jp? / N) 

Hanning: VF(/?)->0.5[l-cos(2*P7/A0] 
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Blackman: W (n) - 0.42 - 0.5 cos(2mi /N) + 0.08 cos(4jih /N) 

Hamming window performed poorly while Blackman window was the best. Figure 26 shows 
the DFT of the previous sampled data after windowing by Blackman window. Each test in this 
report was pre-flltered by Blackman window. Note that the signal is observable only in data set 
#4, where the precession frequency is the smallest. Note also that the windowing attenuates the 
physical resolution by a factor of two. Therefore, although the bin length is still 0.05 Hz, the 
physical resolution is now 0. 1 Hz. 



Frequency [Hz] 



Frequency [Hz] 



Frequency [Hz] 



Frequency [Hz] 


Figure 26 DFT of sampled data after windowing. 


Segmentation 

A hypothetical way to improve the signal detection is to segment the signal with time record 
commensurate with the signal period and to aggregate the segmentations. This technique should 
work well when the white noise is ideal and there is a sufficient number of segmentation such 
that at each time step the white noise will average out. A segmented record should be of an even 
number (since the signal period is 2 sec), but because of our 20-sec total time record, only 
segment records of 2, 5 and 10 seconds are acceptable. Therefore the number of segmentations 
can be only 2, 5 or 10. 


56 






SAO Annual Report #2, NASA Grant NAG3-2881 


x 10 u 



0 0.5 1 1.5 2 



2 3 4 5 

Frequency [Hz] 



12 3 4 

Frequency [Hz] 


Figure 27 Aggregated 10 segmentations for the signal in a white noise. 

Figure 27 shows the segmentation results in the time and in the frequency domains for 10 
segmentations. Since the noise is not pure white and the time record is too short, in addition to 
the signal that adds up, some components of the “white” noise also adds up. Only Monte Carlo 
simulations gave the desired benefit: signal adds up while the uncorrected noise cancelled. 
Another problem is that the resolution deteriorates with increasing number of segments. Since 
we have only few bins separating the signal from the colored noise, reducing the resolution adds 
the noise (colored or white) to the bin related to the signal. Given our current time record and 
sampling rate, it is not recommended to have more than two segments (and probably not to use it 
for a single experiment). In addition, some colored noise may also add up. But only certain 
frequencies, greater than 0.5 Hz and commensurate with the segmentation would add up. Thus, 
integer frequencies from 1 Hz to the Nyquist frequency may add up. 

Amplitude cuts in the frequency domain 

Another possible method is to apply a DFT to the windowed data and then cut or eliminate 
the dominant colored noise in the frequency domain. The advantage of this method is that it 
preserves the signal while reducing the colored noise as needed. Figure 28 shows the data with 
threshold cuts where a white noise was added. The domain of the cuts depends on the particular 
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data. Utilizing a threshold with a magnitude as high as the expected signal can improve the 
numerical convergence of dynamical filters and improve the precision of least square algorithms. 
However, we need to be careful when a white noise is present. The threshold cuts conserves the 
white noise as long as the threshold is above the white noise. But a complete cut eliminates 
certain frequency domain from the while noise, thus destroys the purity of whiteness. 


x 10 16 



xio ' 6 



xio’ 6 x10’ 6 




Figure 28 Windowed data with threshold cuts in the frequency domain. 
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Filtering of colored noise 
Dynamic digital filtering 

The motivation for using digital filters is twofold: To eliminate certain noise frequencies, and 
to resolve the signal. We examined a notch filter and a resonator. Both filters are HR type that 
filter the signal via a dynamical process. The main disadvantage of dynamic filters is the 
transient time for convergence. Therefore, these filters are more adequate for long data spans. 
Nevertheless, we examined the possible advantages of these filters for our data. 

The Notch filter has zeros at the required notch frequency (co denotes the digital frequency). 
The transfer function of the second-order notch filter is: 

HU)-b »-2<*»«y-'+z- 2 

l-2bcosa)nz + ( 2 b-\)z 

U v ' ' 

b — where Act) is the digital frequency width around (O n . 

l+tan(A©/ 2 ) H J o 

The resonator has poles at the required signal frequency. The transfer function of the second- 
order resonator is: 


H(z) - (1 - b) 3- 

l-2&cos<»o£ * + (2&-l)z z 

b - — J — where tan ( ^-) and Gg is the attenuation at a> ■ a>o ± Aco/2. 

1 + /8 -y jl — Gg 2 

In order to make A <o as narrow as possible the radius of the poles (resonator) or zeros (Notch 
filter) must approach the unit circle. However, the filter time constant increases as the radius 
approaches one. Thus, there is a tradeoff between accuracy and convergence speed. The number 
of samples (time constant) to converge to e from the steady state given radius R is 
rj- ln(e)/ln(/?). The relation between the radius and the frequency band is approximately: 
Aa> ■ 2(1- R). The digital frequency for our signal is a>o - 0.0 \n (rad/sample) and the digital 
frequency range, for a frequency band of 0.1 Hz (0.05 Hz from each side), is 
Aa>o ■ 0.002*r (rad/sample). Therefore, the number of samples to converge with c-0.1 is 730 
and the number of samples to converge with e - 0.01 is 1460. Testing a notch filter with our data 
showed three disadvantages. First, some noise shifts toward the signal frequency, especially 
when there is not enough separation between signal and colored noise (data 1 thru 3). Second, 
the transient time is about half of the record time, as explained above. In order to extract the 
signal we need to run the resonator for the rest of the data (the tail) or to apply another signal 
extraction method. The tail is marginally sufficient for another dynamical filter and also contains 
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only about half of the original physical resolution. Third, a notch filter is perfect for removing an 
isolated harmonic noise (for instance, removal of the AC interference from ACG signal) but is 
less efficient for removing a non-stationary noise (data 2 thru 4). 


We tested the resonator as follows. We first pre-filtered the data using amplitude cut filter 
with amplitude threshold of 10 15 m (about three times larger than the expected signal). With 
a>o ■ 0.0 Ur we estimated the amplitude at © 0 using a frequency band of 0. 1 Hz. 






Figure 29 The output of a resonator applied to the raw data. 

Figure 29, that shows the filtered data without pre-processing the cuts, demonstrates the 
necessity of the amplitude cut. The amplitude at 0.5 Hz is an alias of the colored noise at lower 
frequency. We obviously can not reveal the signal with the original colored noise. Note that a 
possible peak at 0.5 Hz is due to the choice of the peaking frequency and not to the real signal. 

Figures 30 and 3 1 show the results from the application of the resonator to the pre-filtered 
noise. The noise in Fig. 30 is colored while the noise in Fig. 31 is a combination of colored and 
white (o= lx 10~ 15 ). 
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Frequency (Hz) Frequency (Hz) 


Figure 30 The output of a resonator applied to a filtered data (no white noise). 
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Figure 31 Output of a resonator applied to a filtered noisy data. 

Figures 32 and 33 show validation tests in which the signal is absent. The noise in Fig. 32 is 
colored while the noise in Fig. 33 is white. It is obvious that the resonator picks some low 
frequency colored noise when the separation of signal/noise is poor. In addition, it picks some 
white noise. Since generally speaking there is no pure white noise, the component of the white 
noise that the resonator picks is always unpredictable. 
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Frequency (Hz) Frequency (Hz) 

Figure 33 The output of a resonator for a white noise (no signal). 

Least square filtering 

Since the signal frequency as well as some components of the colored noise frequencies are 
known, we try to fit these harmonics. The model should include as many known frequencies as 
possible. For optimal fit, the unmodeled noise should be white. In any case, there are some 
unmodeled frequencies, so the estimator is sub-optimal in the best case. One may attempt to 
estimate the colored noise amplitude and then subtract the estimation from the data. This process 
would not work well unless all the modeled colored noise consists of pure harmonics and all 
frequencies are modeled. We modeled the estimator as: X - //0 + n c + n w where X is an N - array 
sampled data, 0 is an 2n + 3array of unknown amplitudes, and n c ,n w are the unmodeled 
colored noise and a white Gaussian noise respectively. There are n modeled frequencies and an 
additional free frequency. The free frequency is close to the signal frequency. The reason will be 
explained later. The corresponding iVx(2n+3) model matrix is: 
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1 cos<»]fi sinajjf] ••• sina)„fj cos(a>„ +fia))f! sin(©„ + 6a>)/] 

1 cos®i ?2 sin®]? 2 ••• sin ®„?2 cos(©„ + 0 a >)?2 sin(a>„ + &»)/ 2 


[l cos a>\t N sin a)\t N ••• sin <o n t N cos (a> n +doo)t N sin(a>„ +d<u)t N 


where (Oj - 2n fj . 


Because of is the relative closeness of some frequencies, we applied a Singular Value 

T 

Decomposition to compute the pseudo-inverse. That is, H - USV where U is N x Vorthogonal 
matrix, V is (2«+ 3) x(2n + 3) orthogonal matrix, and S is Nx(2n+ 3) diagonal matrix that 
consists of the singular values. The estimate amplitude array is: 

0 - V(S T S)~ l S T U T with a Covariance matrix of : 8 - o 2 (H T H) -i . 


We tested the Least Square algorithm as follows. After pre-filtering (as done in the previous 
section) we applied LS with four frequencies (in addition to a bias): two dominants colored noise 
frequencies, the signal frequency and one frequency near the signal. 

We tested the LS with the four data sets and a White Gaussian noise (WGN) with a variance 
o= 10~ 16 (see Fig. 24). We estimated the amplitude (squares) at the signal frequency and the 
amplitudes at a couple of frequencies on each side of the signal (note that the model matrix 
becomes ill defined if a frequency is too close to the signal). The circles show the amplitude of 
the noise without the signal (null signal). The size of the squares is the error estimate. The 
criteria for the existence of a signal should be a visible amplitude peak. Only the signal at data #4 
has a clear peak. Data #3 has some peak while the amplitude estimation of data # 1 and data #2 
does not reveal the signal. If we had a reference noise (circles) we could reveals the signal from 
each data. 
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Figure 34 Amplitude estimation in the vicinity of the signal (with white noise) using LS. 

White noise filtering 
Model 

We assume that we have eliminated the colored noise so the processed data consists of a pure 
white noise and a hypothetical signal. We further assume that the white noise is a WGN with a 
variance cr. 

There are two hypotheses that describe the possible scenarios: 

Noise only: Hg : jc[n]» w[n] 

Noise and a signal: H l : xfn]- s[n]+ 

Our purpose is to detect and estimate the signal. It is not sufficient to estimate the signal 
amplitude because we may pick up a noise with the same frequency as the signal. For example, a 
resonator would show a peak at the desired frequency (for a dense white noise), even for a null 
signal. Therefore, in addition to an amplitude estimator we need to evaluate the confidence of 
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our estimation. Given a probability of false alarm (as small as possible) we would like to 
maximize the probability of detection. 

The main assumption is that the signal frequency and phase are known and the noise is a 
white-Gaussian, therefore: 


(2 no ) 2 a a- o 

1 1 N-l 

P(x\H |)- 2"-/v7T expf 2 2 (4«1- •*[«]) ] 

(2«r ) 2 a a-o 


where x - [x[0] a[ 1] 1]] 7 ^. We tested a unit signal with different SNR. The sampled data 

(dots) is presented in Fig. 35 together with the true signal (solid line). 



time (sec) 


Time (sec) 


Figure 35 Signal with a white noise. 
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Signal detection 

Each sample data x[n ] is considered as a random variable with two possible probability 
distributions: a Gaussian with zero-mean for the null hypothesis, Hq , and a Gaussian centered at 
the expected amplitude for H x . We then assign a threshold y and classify a detection or null- 
detection according to the value of the data with respect to the threshold. The threshold is related 
to the probability of false alarm and to the probability of detection as follows: 

P F A - P(H { ;H 0 ) - P{x > y;H 0 } - jf p(x-,H 0 )dx 

Y 

P D - P(H X ;//, ) - P{x > > - 7 pixMOdx 

Y 

The threshold is a function of the probability of false alarm. Given P FA our aim is to find the 
amplitude that maximizes p d - Note that in general the unknown parameters are amplitude, phase 
and frequency. We used the Generalized Likelihood Ratio Test (GLRT) and decided for H\ 
according the following likelihood criterion: 

pOc 1 A : H x )_ > 

P(x;Hq) 

where A is the Maximum Likelihood Estimation (MLE) of A given in the next section. 

When the amplitude is the only unknown we can express H x as: As[«]+ tv[n), where S is a 
given reference signal. This case is much simpler than the general case (unknown amplitude, 
phase and frequency). 

The model can be written as: x»H>4+ w, where H - [s[0] sfll •••sf/V- 111 7 , therefore the 
hypothesis is now: 

H q : A-0 

H x : A > 0 

It can be shown that we can select H x (for unknown amplitude) if 

N - 1 

Jjc[nWn] N _ { 

fry >y'-a 2 lny + - 2 (4«1) 2 

2 (4«1) 2 M "° 

/ 7-0 

where y' is related to PfA by Pfa “ 2Q(- N /y 7 ). 
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The best P D in the sense of maximum likelihood is: 



2 2 

where the signal energy is defined as: E « A 2 (s[«]) z • 

/?■»() 



Figure 36 Probability of detection Vs energy-to-noise ratio. 


The relation between the signal-to-noise energy ratio (log 10 



) and the probability 


of 


detection, given several probabilities of false alarm is shown in Fig. 36. Our signal is obviously 
detectable with excellent probability for noise variance of few orders of magnitude greater than 
the signal amplitude. However, probability of detection does not tell about the accuracy of the 
signal amplitude estimate (see later). 


Signal estimation 


The amplitude is the only unknown so the amplitude estimate becomes quite simple using the 
model: x m H A + w. The MLE of the amplitude is: 
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N - 1 

24n]s[n] 

A - (H^Hr'H 7 * - -^=2 

N - 1 . 

2 (4«1) 2 

H-0 

We examined the quality of the above estimation for a signal with a variety of white noise. 
The estimation quality depends only on the SNR, so we assign the signal as a unit sinusoid with 
frequency of 0.5 Hz and checked it with different noise variances. Since the 20-s computer 
generated white noise is not a pure white noise, we repeated the test in a Monte Carlo fashion. 
The plots in Fig. 37 show the probability distribution (cumulative frequency distribution) of the 
amplitude estimate for 100 runs. Each plot represents a different SNR. The broken line indicates 
the mean estimate (the true amplitude is one). The probability distribution approach a Gaussian 
as the number of runs increases, as expected from the problem formulation. The amplitude 
estimation is quite good considering that the white noise is greater than the signal. 
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Figure 37 Amplitude cumulative frequency distributions using MLE for different signal-to- 
noise ratios. 
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Summary and Recommendations 

We summarize the advantages and disadvantages of the methods examined in Table 1. 


Table 1 Comparison of the various methods. 


Method 

Advantage 

Disadvantage 

Windowing 

Better observability 

Less resolution 

Segmentation 

Better SNR 

Less resolution 

Amplitude cuts 

Removes colored noise in 
specific band 

Colors the white noise 

Notch Filter 

Remove colored noise 

- Requires transient time 

less resolution 

-The colored noise should be 
in isolated harmonics 

- Does not preserve signal 

Resonator 

Picks the signal directly 

- Requires transient time 

less resolution 

- Picks the white noise 

- No detection confidence 
estimate 

LS Filter 

Estimates the amplitudes of 
the colored noise and the 
signal amplitude 

- Requires a model of the 
noise 

- Only estimate error but no 
detection confidence estimate 

GLRT 

Estimates the amplitude of the 
signal and estimates the 
detection confidence 

Requires a pure white noise 


The method to use for the real experimental data depends on the nature of the data. Pre- 
filtering such as windowing and amplitude cuts are recommended for any type of data. Filters 
such as the resonator and the least square can be used when the separation between the signal 
and the low frequency colored noise is sufficient. Least square filter is preferable for a large 
white noise. However, least square requires good modeling of the noise while the resonator is 
model free. GLRT is the best method when the white noise is large. However, it is optimal for a 
pure WGN while the whiteness of our post-processed data is not perfect. 

The most important future work is to test the performance of the GLRT with a filtered 
colored noise. We executed some preliminary tests of complete cuts that eliminated the noise 
bellow the signal frequency and up to 3 Hz. We tested the GLRT with the inverse Fourier 
transform of the processed data. Except data #1, the estimation is reasonable and very good for 
the case of data #4. Note that although the white noise is far from ideal, the amplitude estimation 
works quite well. Data #4 is a simulated output of the purely rotational detector configuration 
with pivot axis perpendicular to the spin axis. Note also that the colored noise level for the 
configuration with pivot axis parallel to the spin axis is even lower than in data #4. 
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In conclusion, the approach that gave the best results for the extraction of the small EP 
violation signal from a prevalent colored and white noise is as follows: 

Apply Window to data -> Apply Discrete Fourier Transform -> Remove colored noise in 
frequency domain through amplitude cuts ->Apply Inverse Discrete Fourier Transform -> Apply 
Maximum Likelihood Estimation and match filtering. 

The data set #4 (from the purely-rotational detector’s configuration) proved to be the most 
amenable to signal extraction thanks to the good frequency separation between the colored noise 
and the signal, and the low level of dynamics-related noise at the signal frequency. 
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Symbols and Abbreviations 

A = signal amplitude (unknown) 

DFT = Discrete Fourier Transform 

E(x) = expected value of x 

E = signal energy 

FFT = Fast Fourier Transform 

FIR = Finite Impulse Responce 

GLRT - Generalize Likelihood Ratio Test 

H = observation matrix 

H = hypothesis 

I1R = Infinite Impulse Responce 

LS = least square 

Y = threshold 

MLE = Maximum Likelihood Estimation 

N = number of samples 

p(x) = probability density function 
piH^Hj) = probability of H, when Hj is true 
P D = probability of detection 

PfA = probability of false alarm 

Q(x ) = probability that a random variable exceeded x 

in a Gaussian probability distribution 

a 2 = variance 

s = signal 

SNR = Signal to Noise Ratio 

w = noise 

x = data 
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Mechanical Design of Spin and Release Mechanism 

System Operation Description 

The most important requirement on the sensor support and release system is that it releases 
the sensor smoothly during the experiment. In order to do this, the support system must hold the 
sensor against large launch loads, rotate the sensor package up to the desired spin frequency 
(~0.5Hz) and release the sensor without imparting unacceptable lateral dynamics to the package. 
This is accomplished by a design assembled from a number of specialized systems. At the inside, 
the instrument package is held on either side in a 3 point kinematic connection consisting of 3 
ball joints each. The balls are arranged to avoid over-constraint. To ensure that possibly large 
launch loads do not knock the sensor out of position, there is a launch lock for the launch phase 
of the flight. Once aloft, the launch lock will be removed from position. Next there are two (one 
on each side) voice coil/flexure linear stages that allow the sensor to be released smoothly. These 
are each, in turn, mounted on smooth rotary stages that accomplish the sensor spin up. Finally, 
this whole svstem is mounted on a standard hall slide linear stage which pulls this whole system 
out of the way (see Fig. 38). 

Well before release, the sensor package is spun up to 0.5 Hz. Release is accomplished by the 
two equivalent voice coils on either side of the sensor. These pull quickly the ball mounts out of 
position on a set of flexures. The system is designed to have no sliding elements that take part in 
this motion. Since the voice coil/flexure system has a very limited travel, once the sensor is 
rotating freely, the entire sensor support is moved out of the way by the two linear stages that 
support each of the voice coil actuated mounts. The sensor package is now in position to rotate 
without impacting the mass of the supporting structure. 


Mechanical Description - Release Mechanism and Support Structure 
Rotation and Release Assembly 

The instrument package is held in place and rotated between two, three-ball kinematic 
mounts. Each kinematic mount containing the three balls is fastened to a linear voice coil 
actuator. To release the detector package, the actuators are energized in reverse to pull the 
kinematic mounts away by 4 mm each end along the rotation axis. The 4-mm travel will take 0. 1 
sec. The rotating, actuator field assemblies retract the kinematic mounts while the non-rotating 
coil assemblies are attached to the stator housing. See Rotation and Release Assembly, Fig. 39. 
Thin, Titanium flexures are attached to both ends of the field assembly and both ends of the 
internal, rotating actuator housing. A brushless DC motor rotates the instrument package through 
the flexures, voice coil actuators and kinematic mounts. 
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Long Travel Linear Actuators and Support Structures 

Each Rotation and Release Assembly is held in place with side plates attached to a linear 
motion (LM) block (see Fig. 38) mounted to a rail and experiment support channel above the 
instrument package and its clearance envelope. 

Additional linear travel of 10 cm, in about 1 s, immediately after the 4-mm release translation 
is accomplished using an electric thrust cylinder (2 req'd) coupled to the LM block and mounted 
to brackets attached to the experiment support channel. The experiment support channel is 
fastened to mounting brackets welded to the inside of the cryostat ID. 



Figure 38 Spin and release mechanism and its support structure. 
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Figure 39 Spin system and voice coil actuator (first stage of release). 
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Parts specification 

ELECTRICAL THRUST CYLINDER (2 REQ'D): 

IDC MODEL NO. N2-S32 (V) -25-5B (WITH MOTOR BRAKE) 

2.5 / 1 HELICAL GEARS, 5 REV INCH BALL SCREW 
MAX SPEED: 4 IN/S (100 mm/s) 

MAX THRUST: 600 LB (2670 N ) 

LINEAR STROKE: 4 IN. (100 mm) 

MOTOR: 1 .8 deg/step PERMANENT MAGNET HYBRID STEP MOTOR 
POWER (PEAK): 150 W 

POWER INTERMITT ANT DUTY CYCLE: 100 W 
WEIGHT (QTY 2 CYLINDERS): 19.0 LB (8.6 kg) 

LINEAR MOTION GUIDE RAIL AND LM BLOCK (2 RAILS AND 2 BLOCKS REQ'D): 
THK MODEL NO. HRW 50 CA 

4-ROW CIRCULAR ARC CONTACT DESIGN WITH RECIRCULATING BALLS 
DYNAMIC LOAD RATING: 1 1,295 LB (50.2 kg) 

STATIC LOAD RATING: 18,338 LB (81.5 kg) 

MATERIALS (BLOCKS, RAILS AND BALLS): STAINLESS STEEL 
LENGTH OF RAILS: 10.3 IN.(26 cm) 

WEIGHT (TOTAL, RAILS AND BLOCKS): 7 LB (3.2 kg) 

LINEAR VOICE COIL ACTUATOR (2 REQ'D): 

MODIFIED BEI MODEL NO. LA25 - 42 
PEAK FORCE: 60 LB+ (267 N+) 

STROKE REQUIRED: 0.16 IN (4.0 mm) 

STROKE TIME: 0.1 s 
POWER: 375 W 

WEIGHT (QTY 2 ACTUATORS): 4.4 LB (2 kg) 

BRUSHLESS DC MOTOR (2 REQ'D): 

MODIFIED BEI MODEL NO. DIP90 - 15 - 001 Z 
PEAK TORQUE: 3.9 LB-FT (5.3 NM) 

POWER: 28 W 
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EXPECTED RPM @ RELEASE: 30 RPM 
WEIGHT (QTY 2 MOTORS): 1 1 .0 LB (5.0 kg) 

FLEXURES (4 REQ'D): 

MATERIAL: TITANIUM, 6AL4V GRADE 5 
DIMENSIONS: THICKNESS = 0.001 IN (0.025 mm) 

WIDTH = 0.197 IN (5.00 mm) 

LENGTH = 1.181 1 IN (30.00 mm) 

FORCE REQ'D : 50 LB (222 N) 

ROTATION AND RELEASE ASS’Y (HOUSINGS AND BEARINGS): 
MATERIALS: STAINLESS STEEL 

OVERALL DIMENSIONS: 9 3 IN (24 cm ) O.D. X 4.2 IN. (! 1 cm) 
WEIGHT (2 ASSEMBLIES): 20.0 LB (9.1 kg) 

ROTATIONS AND RELEASE ASS'Y SUPPORT STRUCTURE: 
MATERIAL: 6061 -T6 ALUMINUM ALLOY 
WEIGHT (2 ASSEMBLIES): 9.0 LB (4.1 kg) 

THRUST CYCLINDER MOUNTING BRACKET STRUCTURE: 
MATERIAL : 6061 -T6 ALUMINUM ALLOY 
WEIGHT : 21 LB (9.5 kg) 

EXPERIMENT SUPPORT CHANNEL: 

MATERIAL: 6061-T6 ALUMINUM ASSOCIATION CHANNEL 
WEIGHT: 9 LB (4.1 kg). 
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Summary of last-year activity at IFSI laboratory 1 

Construction of new differential accelerometer prototype 

A new differential accelerometer prototype was built by the team at the IFSI laboratory. The 
new prototype consists of proof masses with pivot axes perpendicular to the spin axis (i.e., the 
vertical axis in Fig. 40(a)) and offset from the CM. The proof masses, which are concentric 
hollow cylinders, are attached to the detector’s housing by means of elastic flexures at their 
upper edges (see Fig. 41). Unlike the previous prototype, the capacitive pick-ups are on the 
inner and outer surfaces of the cylinders. 



Figure 40 New differential accelerometer prototype: (a) cut-up view of proof masses and signal 
pickups; and (b) differential accelerometer assembled inside its vacuum chamber. 

This pickup arrangement provides for a better electro-mechanical coupling of the 
accelerometer (i.e., a better conversion of mechanical energy into electrical energy) and an 
increase of the effective mass of the proof masses. Moreover, in the new prototype the CMs of 
the proof masses are coincident within the attainable construction accuracy. 

1 Section contributed by V. Iafolla, S. Nozzoli, D. Lucchesi, F. Santoli of the Institute of Interplanetary Space 
Phyiscs (IFSI) in Rome, Italy. 
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Figure 41 Finite element analysis of the proof mass vibration modes. 

The proof masses in this prototype are both made of aluminum. The internal mass is 0.96 kg 
and the external mass is 1.1 kg. The resonant frequency in a 1-g environment depends on the 
top-up (regular pendulum) or top-down (inverted pendulum) orientation of the detector. The 
resonant frequency for the top-up configuration (that is normally used in the laboratory) is 2.28 
Hz. Table 2 summarizes the key characteristics of the new differential accelerometer. 


Table 2 Key characteristics of differential accelerometer prototype. 


ITEM 

VALUE 

Internal proof mass 

0.96 kg 

External proof mass 

11kg 

Resonant frequency (top-up) 

2.28 Hz 

Pre-amplifier noise temperature 
(commercial amplifier) 

0.75 K 


The development of the new prototype has occurred in parallel to the progress made in the 
analysis of the configuration for the flight-experiment detector. The new prototype adopts the 
roto-translational configuration which has many features in common with the purely rotational 
configuration. On the one hand, the roto-translation does not have the same degree of rejection of 
dynamics-related noise of the purely rotational configuration. On the other hand, the roto- 
translational configuration is amenable to more meaningful testing on ground because the motion 
of its test masses can be excited gravitationally by tilting the detector off the local vertical while 
the purely rotational configuration is practically insensitive to tilt. 
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Noise measurements of the new detector prototype 
Introduction 

This section illustrates the results of the experimental activity carried out in 2004 to evaluate 
the noise level of the new differential accelerometer prototype. 

Measurements have been carried out with four different experimental setups: 

1) The fist setup consists of the measurement electronics coupled to two reference 
capacitances on each channel that simulate the characteristics of the accelerometer. Lock- 
in amplifiers have been used for generating the polarization potential, signal 
demodulation and acquisition. This setup has been used to evaluate the effect of 
atmospheric pressure on the noise level of the electronics. 

2) The previous electronics setup has been coupled to one sensor of the differential 
accelerometer to measure the seismic noise level. Moreover, we evaluated the input 
noise to the lock-in amplifiers. 

3) Amplifiers have been added to the previous setup to amplify the polarization tension and 
increase the signal-to-noise ratio. We carried out measurements with reference 
capacitances on one channel while the other channel was coupled to the accelerometer 
sensor (excited by seismic noise). 

4) In this setup, we added the input pre-amplifiers to the lock-in amplifiers to increase the 
signal-to-noise level and we carried out measurements of the electronics-related noise on 
one channel (with reference capacitances) and seismic noise measurements on the other 
channel. 

The data acquisition was carried out with a sampling frequency of 64 Hz with data strings of 
about 15-minute each. A low-pass filter with a time constant of 10 ms and a cut-off attenuation 
of 1 4 dB/oct was used. 


Summary of experimental results 

Here we show a short summary of results while more details can be obtained from the 
relevant documents written by the IFSI laboratory team. 

Table 3 shows a summary of the output voltage noise levels measured with setup 1 in which 
the measurement electronics is coupled to reference capacitances (simulating the acceleration 
sensors). 
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Table 3 Noise levels measured with setup 1 . 


Channel 

Acquisition 

Pressure 

(mBar) 

Noise at 1 Hz 
Wsqrt(Hz) 

Noise at 10 2 Hz 
V/sqrt(Hz) 

1 

EL 02 

ambient 

2X1 O' 6 

2x10 5 

1 

EL 03 

6.7x1 0' 3 

lxlO' 6 

9x1 O’ 6 

1 

EL 05 

5.9X10 3 

2x1 O' 6 

4x1 O' 6 

2 

EL 02 

ambient 

2x1 O' 6 

2x1 O' 5 

2 

EL 03 

6. 7x1 0' 3 

2x1 O' 6 

2x1 O' 6 

2 

EL 04 

5.9X10 3 

7x1 0' 7 

1X10 6 


In Setup 2 one channel was connected to the internal proof mass of the accelerometer while 
the other channel was connected to the reference capacitances. The two channels are marked by 
the letters I and E which are associated with the internal and external proof masses of the 
differential accelerometer. The measurements with the lowest noise level for this setup are 
shown in Figs. 42(a) and 42(b) that portray the time-domain signal and its frequency spectrum, 
respectively. 

Setup 2 was also used to convert the voltage measurements into acceleration units (g) by 
tilting (by known values) the base plate of the detector prototype. The results of this calibration 
measurements are shown in Fig. 43. 



(cont’d) 
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Finally, we have collected in Fig. 44 a consistent set of results obtained from measurements 
carried out with setup 2 (shown in black), setup 3 (green) and setup 4 (magenta) from the two 
channels. Note that setup 1 was not used to carry out measurements of seismic noise. The 
results are expressed in terms of acceleration power spectral densities ( gl^Hz )■ 

The upper group of lines is associated with channel 1 (connected to the internal proof mass) 
while the lower group of lines is associated with channel 2 (connected to the reference 
capacitances). The upper group of lines represent the seismic noise measured by the internal 
proof mass of the differential accelerometer. These lines are practically alike as it should be 
expected because the seismic noise does not change substantially with time. 

The lower group of lines represent the noise level of the entire electronic measurement chain. 
The noise at low frequencies is attributable to external thermal effects while the noise at higher 
frequencies is associated with the input noise of the lock-in amplifiers. 


g/ 

Hz 1 ' 2 



Frequency (Hz) 


Figure 44 Acceleration spectral densities of a consistent set of measurements taken with setup 2 
(black), 3 (green), and 4 (magenta). The upper curves are seismic noise measurements from the 
channel connected to the internal proof mass of the detector. The lower curves are noise levels 
of the electronic measurement chain from the channel connected to the reference capacitances. 
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The experimental results demonstrate that the electronics-related noise is much lower than 
the seismic noise even with commercial electronics operating at room temperature. 
Consequently, the accuracy of the measurements conducted in the laboratory will be limited by 
the seismic noise. Moreover, the arrangement of setup 4 exhibits the lowest level of electronics- 
related noise of those explored. 


Preliminary measurements of common-mode rejection 
Measurements with one lock-in amplifier 

The measurement (GRT53) was carried out by utilizing the measurement scheme presented 
in Fig. 45. The seismic noise, measured by the internal mass, is sent to one lock-in amplifier and 
then to two channels of the data acquisition system. This measurement is needed to determine 
the upper limit of the common-mode rejection factor achievable with two lock-in amplifiers and 
two sensors. 



Figure 45 Scheme to measure the common-mode rejection of the acquisition system. 

Figure 46 shows the spectral densities of the seismic signal and the rejection obtained at the 
output of the data acquisition system. After comparing the two spectra, we computed a common- 
mode rejection factor in the range 3x10 4 and 3x10 5 in the frequency range 0.008 Hz and the cut- 
off frequency of 1 Hz (see Fig. 47). 
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Figure 46 Spectra of seismic noise acquired by one proof mass and the rejection residual from 
the data acquisition system. 



Figure 47 Common-mode rejection factor associated with data acquisition system. 
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Measurements with two lock-in amplifiers 

With this measurement set (GTR54) we aim to evaluate the contribution to the common- 
mode rejection of the two lock-in amplifiers. The seismic signal measured by the one proof 
mass has been sent to two lock-in amplifiers and then to the data acquisition system as shown in 
Fig. 48. 



V out 
1 


V out 
2 


Figure 48 Arrangement to measure the common-mode rejection of the two lock-in amplifiers. 

The input (seismic) signal vs. time is shown in Fig. 49. The spectra of the signal and of the 
rejection obtained at the output of the measurement chain are shown in Fig. 50. 
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Figure 49 Input signal vs. time. 



Figure 50 Spectra of input signal and output rejection signal. 
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After comparing the two spectra we compute a common-mode rejection factor ranging from 
100 to 3000 in the frequency band 0.007 Hz and the cut-off frequency of 1 Hz (see Fig. 51). 



Figure 51 Contribution to common-mode rejection of the two lock-in amplifiers. 

Calibration in acceleration units and rejection factor 

A new calibration curve, to be associated with this measurement chain, was derived (see Fig. 
52) to obtain a calibration factor for the internal proof mass of a 2. 149 10 3 Volt/g. 



Figure 52 Calibration curve Volts - g for the internal proof mass. 
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By adopting the calibration factor computed previously, the spectra of the input signal and 
rejection output can be expressed in terms of acceleration units. The spectra, depicted in Fig. 53, 
show that the lowest rejection value is 10 10 g/sqrt(Hz) for frequencies below the cut-off 
frequency of 1 Hz. This value must be considered the limiting value for this particular 
measurement chain. 



Frequency (Hz) 

Figure 53 Spectra of input signal and rejection output in terms of acceleration units. 


As previously stated, these are preliminary measurements of the common-mode rejection 
factor in which the limit imposed by the noise of the electronics measurement chain has not yet 
been approached. 
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Conclusions 

The analysis of the dynamics in free fall of various detector configurations in the presence of 
release errors and proof masses imperfections has yielded the design of a suitable detector 
configuration. This configuration provides a very low level of dynamics-related noise at the 
signal frequency that meets the deired detector’s threshold sensitivity goal. The configuration 
consists of proof masses that have purely rotational degrees of freedom. The internal proof mass 
is the “EP-violation-sensitive” mass that, thanks to its design made of two halves of different 
materials, senses an EP violation as a torque about its pivot axis. The second mass, made of one 
of the materials of the other mass, is the “dynamical-reference” mass. The role of this proof 
mass is to attenuate (and ideally remove), through output differencing, the dynamics-related 
noise from the detector’s output. In this design, an EP violation will produce a torque that is 
comparatively strong with respect to gravity gradient torques. 

We have demonstrated through extensive simulations that this detector design is capable of 
detecting a very small differential signal and resolving that signal out of the dynamics-related 
noise and gravity gradients in the presence of conservative values of the construction 
imperfections of the proof masses and errors at release. This conclusion implies that, with this 
new design, the dynamics-related noise is no longer the limiting factor of the detector's 
threshold sensitivity. The limiting factor becomes the intrinsic noise of the detector (e.g., the 
Brownian noise of the proof masses). 

In addition, we have devised a data analysis technique that enables the extraction of a small 
signal from the colored noise (associated with the detector’s dynamics and gravity gradients) and 
white noise (associated with the detector’s intrinsic noise). Moreover, we have estimated the 
probability of detection of a monochromatic signal of known frequency and phase and the likely 
error in the amplitude estimate as a function of the signal-to-noise ratio. 

Our engineering department has developed a design of the mechanism to spin and release the 
instrument package inside the cryostat. The mechanism will be capable of operating at cryostatic 
temperature. The mechanism spins up the detector and then releases it through a two stage 
process: (1) a very quick short retraction of the two kinematic mounts that support the instrument 
package; and (2) a relatively fast and longer retraction of the first stage to clear the mechanism 
out of the instrument package motion. 

Our non-US partners have built a second differential detector prototype that incorporates 
many of the features planned for the flight-experiment detector. They also carried out 
preliminary laboratory measurements of instrument noise and common-mode rejection factor of 
the new detector protoype. 
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Planned Future Work 


Tasks to be carried out during the next year research activity are as follows: 

Design and construction of an engineering model of the spin and release mechanism for room 
temperature testing. The engineering model will have all the functionality of the (cryogenic) 
mechanism for the flight experiments and provisions for unloading the weight of the instrument 
package so that it can be used for laboratory tests that mimic the zero-g release dynamics. 

Compute higher-order gravitational harmonics (coupled to the capsule mass) of the internal 
proof mass that is made of two halves of different materials, in order to verify that the new 
design does not violate the noise level requirement of the detector. 


Apply the newly-developed noise extraction techniques to data strings obtained from the 
laboratory noise measurements of the new differential accelerometer prototype. A signal with a 
signal-to-noise ratio (SNR) comparable to what is expected in the flight tests will be added to the 
noisy data to test the signal-extraction strategy on experimental data. 


Derive a preliminary cost estimate for the design and construction of the balloon-released 
facility and the recurring operating costs for the experiments. 


Provide inputs to our non-US partners for a detailed design of the flight detector based on the 
results of the detector’s dynamics response analysis conducted thus far. 


93 



SAO Annual Report #2, NASA Grant NAG3-288 1 


Papers published or submitted during last year grant activity 

Ashenberg, J. and E.C. Lorenzini, “Analytical formulation of a complex mutual gravitational 
field.” Classical and Quantum Gravity. Vol. 21, No. 8, 2089-2100, 2004. 

Iafolla, V., D.M. Lucchesi, S. Nozzoli, M. Ravenna, F. Santoli, 1. 1. Shapiro, E.C. Lorenzini, 
M.L. Cosmo, C. Bombardelli' J. Ashenberg, P.N. Cheimets, and S. Glashow, “General 
Relativity Accuracy Test (GReAT): new configuration for the differential accelerometer.” 
Presented at the COSPAR-2004 Congress, Paris, France, July 2004. Also submitted to 
Advances in Space Research . 

Lorenzini, E.C., 1. 1. Shapiro, M.L. Cosmo, P.N. Cheimets, C. Bombardelli, V. Iafolla, D.M. 
Lucchesi, S. Nozzoli, M. Ravenna, F. Santoli, and S. Glashow, “Development of the 
TEPEE/GReAT experiment to test the Equivalence Principle in Free Fall from a balloon.” 
Proceedings of Testing the equivalence principle in Space and on Ground, Pescara, Italy, 
20-23 September 2004. 

Iafolla, V., D.M. Lucchesi, S. Nozzoli, M. Ravenna, F. Santoli, I.I. Shapiro, E.C. Lorenzini, M.L. 
Cosmo, J. Ashenberg, P.N. Cheimets, and S. Glashow, “The Weak Equivalence Principle 
(WEP) and the General Relativity Accuracy Test (GReAT) with an Einstein Elevator.” 
Proceedings of 16 th SIGRAV Conference on General Relativity and Gravitational Physics, 
Vietri, Italy, 13-16 September 2004. 


94 


SAO Annual Report #2, NASA Grant NAG3-2881 


References 


' Shapiro, I.I., et al. “Test of the Equivalence Principle in an Einstein Elevator.” Smithsonian Astrophysical 
Observatory, Final Report on NASA Grant NAS8-1780, April 2004. 

11 Baessler, S., B. Heckel, E. Adelberger, J. Gundlack, U. Schmidt, and E. Swanson, Phys Rev Lett., Vol. 83, 1 
Nov. 1999. 

“ Williams, J.G., X.X. Newhall and J.O. Dickey, Relativity parameters determined from lunar laser ranging, Phys. 
Rev. £), 53, p. 6730, 1996. 

,v Damour T, F. Piazza, and G. Veneziano, “Runaway dilaton and equivalence Principle violations.” Physical 
Review D, Vol. 66, 081601-1, 2002. 

v Damour T, F. Piazza, and G. Veneziano, “Violations of the Equivalence principle in a dilaton-runaway 
scenario.” Physical Review Letters Vol. 89, 046007-1, 2002. 

V1 Lorenzini, E.C., I.I. Shapiro, F. Fuligni, V. lafolla, M.L. Cosmo, M.D. Grossi, P.N. Cheimets and J.B. Zielinski, 
“Test of the Weak-Equivalence Principle in an Einstein Elevator.” II Nuovo Cimento, Vol. 109B, No. 11, 1994. 
vii Parzianello, G., E.C. Lorenzini, I.I. Shapiro, M.L. Cosmo, V. lafolla and S. Nozzoli, “Analysis of 
Effects of Initial Conditions and Physical Parameters in a Differential Accelerometer Modeled as a 
Multibody System.” Presented at the 54 ,h International Astronautical Federation (IAF) Congress, Oct, 

n r* 

z,wj , uicincii, wciiiiaiiy. 

v,u Shapiro, I.I., et al. “Test of the Equivalence Principle in an Einstein Elevator.” Smithsonian Astrophysical 
Observatory, Annual Report #1 on NASA Grant NAS3-2881, April 2004. 

,x lafolla, V., A. Mandiello and S. Nozzoli, “The High Sensitivity Italian Spring Accelerometer (ISA) for 
Fundamental Physics in Space” Advances in Space Research . Vol. 25, No. 6, 1241-1244, 2000. 


95 



